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ABSTRACT 

Luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) dominate the star formation 
rate budget of the universe at z > 1, yet no local measurements of their heavy element abundances 
exist. We measure nuclear or near-nuclear oxygen abundances in a sample of 100 star-forming LIRGs 
and ULIRGs using new, previously published, and archival spectroscopy of strong emission lines 
(including [O II]AA3727, 3729) in galaxies with redshifts (z) ~ 0.1. When compared to local emission- 
line galaxies of similar luminosity and mass (using the near-infrared luminosity-mctallicity and mass- 
metallicity relations), we find that LIRGs and ULIRGs are under-abundant by a factor of two on 
average. As a corollary, LIRGs and ULIRGs also have smaller effective yields. We conclude that the 
observed under-abundance results from the combination of a decrease of abundance with increasing 
radius in the progenitor galaxies and strong, interaction- or merger-induced gas inflow into the galaxy 
nucleus. This conclusion demonstrates that local abundance scaling relations are not universal, a 
fact that must be accounted for when interpreting abundances earlier in the universe's history when 
merger-induced star formation was the dominant mode. We use our local sample to compare to high- 
redshift samples and assess abundance evolution in LIRGs and ULIRGs. We find that abundances in 
these systems increased by ^^0.2 dex from z ^ 0.6 to z ~ 0.1. Evolution from z ~ 2 submillimeter 
galaxies to 2: ~ 0.1 ULIRGs also appears to be present, though uncertainty due to spectroscopic 
limitations is large. 



Subject headings: galaxies: abundances — galaxies: evolution — galaxies: interactions 
ISM — galaxies: kinematics and dynamics — infrared: galaxies 



galaxies: 



1. INTRODUCTION 

Mid-infrared and submillimeter observations show that 
luminous and ultraluminous infrared galaxies (LIRGs 
and ULIRGs)^ h ost most of the star formation in the 
universe at z > 1 (Le Floc'h et al. 2005: Chapman et al 



20051 iDaddi et al.H2005l: IWang et al.H2006l : ICaputi et al' 



quasar is left ([Sanders et al.l Il988f l. This hypothe- 
sis is under scrutiny; comparison to the fundamental 
plane of ellipticals shows that many ULIRGs are evolv- 
ing into moderate-mass elhpticals (iGenzel et al.l 120011: 
Veilleux et all 120021 : iTacconi et all 120021 : iDasvra et all 



20071 ). Understanding local examples of these sources is 



thus a window to star formation and accretion onto su- 
permassive black holes at the epoch of highest star for- 
mati on rate and a,ctive galactic nucleus (AGN) density 
(e.g.. lMadau et al.llT998l : ISchmidt et al.| [r995l. 

A great deal is known about nearby ULIRGs; for a 
recent review see iLonsdale et all (l2006h . All ULIRGs 
possess strong starbursts, and many also host optically- 
visible AGN (Veilleux ct al. 1999). The starburst s and 
AGN are on averag e heavily obscured (Genz el et al.l 
119981 iHao et al.l[2007fl. and o p tically-invisible AG N may 
be present (|Lutz et al.lll999l : lArmus et al.ll2007f ). The 
prevalence of optically-visible AGN increases with in- 
creasing infrared luminosity (Veilleux ct al. 1999) and 
as the merger progresses (Veill eux et al.l 120021 ). It is 
thus hypothesized that ULIRGs play a role in the 
evolution of quasars. When obscuring dust is re- 
moved from a buried AGN (by st arburst- or AGN- 
driven outflows; see, e.g.. lRupke et al. 2005a, c). a bright 



2006alll J)T similar to but slightly less massive than the 
hosts of optically bright quasars (Dasyra et al. 200^. 

LIRGs host less intense starbursts than ULIRGs, and 
the frequency of occurr ence of optically-vis ible AGN in 
LIRGs is much smaller (IVeilleux et al.lll995[). Photomet - 
ric studies divide LIRGs into two groups llshidall2004) . 
The most luminous are early in the merger sequence of 
two roughly equal-mas s galaxies, and thus may be th e 
progenitors of ULIRGs (|Arribas et al.ll2004l : IIshidall2004D . 
Other LIRGs are unequal-mass mergers or isolated disk 
galaxies which may or may not be experiencing an inter- 
action (Ishida 2004). Kinematic studies confirm that the 
pair mass ratios in at least some LIRG interactions are 
of order 1 — 3 (jRothberg fc JosephI 120061 : iDasyra et al.l 
2006a). 

The interstellar medium in LIRGs and ULIRGs is 
in a ki nematically extreme s tate, d ominated by inflows 
'e.g., iBarnes fc HernquTstI 119961: iMihos fc Hernaui"st 



19961) outflows (iHeckman et al.l |2000t iRupke et al 
2002 . I2005al lbll3: iMartinI 120051 120061: see [Veilleux et a. 
20051 for a recent review), and turbulent motions 
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1 LIRGs are defined^byiO^^ < Lir/L© < lO^^ and ULIRGs by 
10^^ < Z/ir/ Lq < lO'^^ , where Ljr is the 'total' infrared luminosity 
from 8 - 1000 ^lm. 



pownes fc Solomon! Il998f ). These gas motions have 
the power to significantly alter the chemical states 
of the progenitor galaxies (e.g., lEdrnun di II990I : 
IKoppen fc Edmundsl li999l: IDalcantonI l2007f) . Ongoing, 
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intense star formation in LIRGs and ULIRGs is also 
producing and redistributing heavy metals at a prodi- 
gious rate. In this paper we describe the first compre- 
hensive study of the oxygen abundances of local LIRGs 
and ULIRGs (§il2] -TB]). This will allow us to assess 
the effects of these processes on the gas-phase abun- 
dances of infrared-selected, interacting galaxies. Studies 
of optically-selected mergers su ggest that g as motions do 
alter nuclear abundances (Kewl ev et al.l f2006a). Here we 
present evidence that this is true also of strong mergers 
with high star formation rates. 

In order to understand the chemical evolution of LIRGs 
and ULIRGs, it is crucial to compare these sources to 
weakly- or non-interacting galaxies with modest star for- 
mation, which represent the progenitors of LIRGs and 
ULIRGs. To this end, we compare LIRGs and ULIRGs to 
published luminosity-metallicity, mass-metallicity, and 
mass-effective yield relations (§S]-t6]). We discuss these 
results in ^J7l 

There exist a few measurements of oxygen abundances 
in hig h-redshift infrared-luminous galaxies. Liang et al.l 
(|2004l) measure the abundances of ^20 z = 0.4 — 0.9 
LIRGs selected at 15 /J,m by the Infrared Space Obser- 
vatory (ISO). Abundance measurements also exist for a 
handful of z ~ 2 submillimeter-selected galaxies that 
have total infrared l uminosities greater than or equal 
to those of ULIRGs (|Tecza et al.ll2004l : ISwinbank etaH 
[200l . Finally, we describe in this paper a handful of 
new moderate-redshift measurements. 

While these high-z measurements are valuable on their 
own, they are also sparse and, in many cases, uncer- 
tain. The more robust measurements of local LIRGs and 
ULIRGs that we report here facilitate two comparisons: 
we can now compare the enrichment histories of infrared- 
luminous and infrared-faint galaxies in the local universe 
(§i) l4.2ll473l and[5]); and we can look for evidence of chem- 
ical evolution among infrared-luminous populations as a 
function of redshift ( ij4.4p . 

We summarize our work and discuss its consequences 
in ga 

Throughout the paper, we use the notation (O/H) to 
refer to the ratio of the number densities of O and H 
atoms in the ISM. The variable Z refers instead to the 
mass fraction of O relative to the total mass of gas. These 
two variables are related by a constant: Z = 16/C(0/II), 
where C ~ 1.4 is the ratio of total to H gas masses. 
Where appropriate, we use (O/H) and Z interchange- 
ably. For the solar oxygen abun dance, we use the recent 
value from lAsplund et all (|2004f) : 12-hlog(0/H)Q = 8.66. 
For all calculations, we assume the standard cosmology 
of Ho = 75 km s'^ Mpc-\ flm = 0.3, and flA = 0.7. 

2. SAMPLE SELECTION 

Abundance diagnostics of star-forming galaxies rely 
mostly on emission lines. The emission lines of many 
LIRGs and ULIRGs, however, may inclu de contributions 
from an AGN and /or strong shocks ( Kim et aTl Il998t 
I Veilleux et al.lfT999f l. Abundance diagnostics are gener- 
ally calibrated using galaxies with modest star formation, 
whose emission line fluxes do not contain strong contribu- 
tions from cither of these ionization mechanisms. Thus, 
choosing galaxies whose lines are starburst-dominated is 
important for computing accurate abundances. 

This decision is complicated by the multiple options 



for defining a star-forming galaxy in the phase space 
of optical emission line fl ux ratios. For instance, usiii g 
the classic diagnostics of IVeilleux fc Osterbrockl (|1987D . 
60% of LIRGs and one-third of all ULIRGs ar e 'H II- 
region- like' galaxies (jVeilleux et al."l995l Il999l ). More 
recent work updates this empirical classification scheme 
using ^10 ^ galaxies from the Sloan Digital Sky Sur - 
vey fSDSS: [Kauffmann et al.l[200llKewlev et al.ll2006bD . 
Many LIR Gs and ULIRGs classifie d as H II galaxies using 
the V eilleux fc Osterbrockl ([1987') scheme lie away from 
the bulk of local star-forming galaxies; instead, they lie 
in the region of the [O ni]A5007/H/3 vs. [N II]A6583/Ha 
flux ratio diagram that is be tween the outer b oundary of 
the locus of SDSS galaxies (jKauffmann et all [20031 and 
the line delineating the maximum li ne ratios a chievable 
by starbursts, according to theory ( Kewlev et al . 2001). 

In ij3.31 we discuss in detail the abundance uncertain- 
ties that arise from such physical effects. For now, we 
adopt loose initial selection criteria. We include in our 
sample all galaxies classified as H II galaxies or low ion- 
ization nuclear emission-line regions (LINERs) under any 
scheme. We also require that line flux uncertainties be 
smaller than 50%; almost all fluxes are far more certain 
than this, but a few fall near the limit. In Figure [1] 
we place the galaxies in the current sample on several 
line-ratio diagrams, with various classification schemes 
over-plotted. 

Our ULIRGs are primarily from the 1 Jy sample, 
which is a complete, flux-limited, northern-hemisphere 
sample drawn from the Infrared Astronomical Satellite 
(IRAS) database (Kim & Sanders 1998). We require 
that our galaxies have measured [O II]AA3727, 3729 
fluxes. T h e moderately -high-resolution spectra of 
iRuDke et all (|2(M l2005bh . from Keck or the MMT, 
have broad enough wavelength coverage for this pur- 
pose. We supplement this data set with a handful 
of spectra of 1 J y objects f rom the Fifth Data Re- 
lease of the SDSS (lYork et al.ll2000t I Strauss et al.l [200l 
lAdelman-McCarthv et all 120071 )" To improve our statis- 
tics, we also include a few ULIRGs with published [O II] 
fluxes from the Revised Bright Galaxy Sample (RBGS), 
the Warm Galaxy Survey (WGS ), and the 2 Jy survey 
()Kim et all 119951: iWu et al.i[l99l . The final local sam- 
ple of 31 galaxy nuclei has (z) = 0.14, with a maximum 
redshift of 0.27. However, in addition we include mea- 
surements for a few g alaxies with z = 0.4 — 0.5 from the 
FIRST-FSG catalog (jStanford et all [20001 : ImTpke et all 
I2005b3 for the purpose of assessing redshift evolution 

(SH). 

We have selected our LIRGs primarily from the 
Revised Bright Galaxy Sample, whi ch is also a 
comp lete, flux- limited IRAS sample (iSande rs et al.l 
|2003[) . Our RBGS data c ome from a variety of 
spectroscopic su r veys: (1) iKim et all (|1995) ; (2) 
Liu &: Kennicup (^1991; (3) IWu et al.l ^M ): (4) 
Rupke et al.l ~([2002. 20(3 53); (5) the Fourth Data Re- 
lea se of the SDS S ( Adclman- Mc Carthy et"al] |2006^: and 
(6') lMoustakas fc Kennicutt (200^ nuclear spectra only). 
Again, to improve our statistics, we include published 
fluxes of galaxies from the Warm Galaxy Survey (^ WGS) 
and 2 Jy sample (|Kim et al.lll995HWu et al.|[T99l . The 
sample of 65 galaxy nuclei has (z) = 0.04. (We also have 
measurements for one LIRG with z = 0.48, which we use 
in SH) 
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Fig. 1. — Emission-line ratio diagrams, with fluxes corrected for extinction. Blue stars are luminous infrared galaxies and red circles 
are ultraluminous infrared galaxies. The solid lines separate star-forming galaxies from LINERs fVeilleux & Osterbrock 1987); the dotted 
lines denote the phase space li mits of the [Kewlev et al.i t2001j ) starburst models in all plots except the bottom right, where they separate 
star bursts, L INERs, and AGN l|Kewlev et ani2006bl) : and the dashed line denotes the outer boundary of the bulk of star- forming galaxies 
in the SDSS I IKauffmann et al.ll2003fK 



The selected galaxies are representative of the local 
{z < 0.2), star-forming, infrared-luminous galaxy pop- 
ulation. The 1 Jy sample, 2 Jy sample, RBGS, WGS 
are complete samples. The particular galaxies which 
appear in this paper are effectively a random selection 
from these samples, since they are culled from their par- 
ent spectroscopic studies only on the basis of spectral 
type and signal-to-noise ratio. The parent studies from 
which our spectroscopic data were collated impose vari- 
ous constraints on sky location, emission-line sensitivity, 
and (to a minor degree) redshift. However, these have 
no effect on the physical properties of the galaxies cho- 
sen, as shown by comparisons among F igure [1] and simi- 
lar figures from the parent studies (e.g.. lKim et al.iri995l : 
iVeilleux et aHflQQflh . 

Some of the galaxies in our sample are multiple-nucleus 



systems^. For the 1 Jy ULIRGs, the nuclei are charac- 
teriz ed using sensitiv e optical and near-in frared imag- 
ing (|Kim et all |200^ : IVeilleux et a"ni2002f) . For other 
sources, we specify nuclei using unique designations from 
the NASA/IPAC Extragalactic Database (NED; see Ta- 
ble [1]). We treat each nucleus separately in our analysis, 
since the apertures for our spectra are nuclear or near- 
nuclear ( ij3.3p . Accordingly, for each multiple- nucleus 
system we divide the total system infrared luminosity be- 
tween component n uclei based on resolved IRAS imaging 
(jSurace et al.l[2004h or the near- infrared luminosities of 
the component galaxies f ij4.ip . The fraction of the to- 
tal infrared luminosity of the system belonging to each 
nucleus (or equivalently fractional near-infrared luminos- 

^ For some multiple-nucleus systems, only one nucleus enters our 
sample. This is due to a lack of data or because the other nucleus 
has a Seyfert optical spectral type. 
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ity) for those systems where the NIR luminosity is used 
to divide iiR is hsted in Table [H Seven ULIRG nu- 
clei descend into the LIRG category due to their nuclear 
luminosity. This reclassification does not impact our re- 
sults. 

In total, there are 100 galaxies or nuclei in our sam- 
ple. Table [1] lists the basic properties of each galaxy or 
nucleus. 

The strong emission lines in the iRupke et all (|2002l 
|2005b) and SDSS flux-calibrated spectra were measured 
using the IRAF task SPLOT. The average measurement 
uncertainty in the brightest hues (e.g., [O II] or [N II]) is 
5% or less. For weak or noisy lines ([O I], [O III], or [S II] 
in a few cases) or those affected by a continuum that 
has strong stellar absorption, the uncertainty rises to ~ 
20 — 30%. Table [5] lists the line fluxes newly measured 
for this study. 

Our Keck, MMT, and SDSS spectra are typically of 
high enough resolution for us to fit simultaneously a 
Voigt absorption and Gaussian emission component to 
H/3. In the few cases where this was not possible, we used 
lower order Balmer lines to estimate the expected absorp- 
tion in H/3. We corrected the Ha emission line in these 
galaxies for stellar absorption by extrapolating from the 
absorption equivalent widths of lower order Balmer lines, 
using the relative value s for different lines predicted by 
the oscillator strengths ljMenzellll969t ): 



VFeg(Ha) 



(1) 



In this equation, H^r is the Balmer transition from which 
the equivalent width of Ha is to be calculated; H^ repre- 
sents the Balmer transition with upper principal quan- 
tum number i (e.g., H4 = H/3); a nd /(H^) is its oscil- 
lator strength. iPatris et al.l ()2003f ) tabulate the coeffi- 
cient values that we used. Table [2] lists the (absorption- 
uncorrected) Ha emission line and H/3 absorption line 
equivalent widths for newly measured data. 

The emission lines were corrected for extinction using 
the Balmer decre ment. We assume an intr insic Ha/H/3 
flux ratio of 2.86 (|Hummer fc Store"vlll987| ). an effective 
foreground scr e en, an d the starburst attenuation curve of 
iCalzetti et all (|2000l ). Table [D lists the E{B-V) values 
for newly measured data. 

The data quality was checked through cross-correlation 
of objects common to different data sets. The results 
showed remarkable consistency, given the number of ref- 
erences from which the data were drawn. The typical dis- 
crepancy between two measurements of the same galaxy 
nucleus is of the order ^0.1 dex. Where multiple data 
were available for a given source, we chose the spectrum 
with highest spectral resolution and signal-to-noise. Ex- 
ceptions are sources not drawn from the 1 Jy sample or 
RBGS that have both published and SDSS spectra; for 
consistency, we chose the published data even though the 
SDSS data may be of higher spectral resolution or sen- 
sitivity. We decided to use the SDSS data only to sup- 
plement the number of spectra available from the 1 Jy 
sample and RBGS (or for data quality checking) in order 
to keep the sample selection fairly clean and the size of 
the sample manageable. (I.e., there are a large number 
of infrared-lumi nous galaxies in the SDSS that we did 
not include; e.g.. iPasquah et aLll2005l ) 



3. ABUNDANCES 
3.1. Comparison of Diagnostics 

Numerous strong- line abundance diagnostics exist, 
each relying on different combinations of the strongest 
emission lines measurable in optical spectra. Each diag- 
nostic in turn has different absolute calibrations, based 
on photoionization models, electron temperature (Tg) 
measurements (i.e., weak-line diagnostics), or a combi- 
nation of the two. Different diagnostics, or different cali- 
brations of the same diagnostic, can give vastly different 
abundances for the same galaxy or group of galaxies. For 
instance, models calibrated on measurements of electron 
temperatures of H II regions in nearby galaxies differ 
from theoretical calibr ations by factors of a few (see, e.g., 
iKennicutt et al.ll2003l and references therein). 

To help the reader understand the uncertainty that at- 
taches to the choice of a particular diagnostic/calibration 
combination, we have investigated several d ifferent op- 
tions. One of these fthe lTremonti et "ani200 4 calibration') 
we employ in §§|3]— [Elto compare our data with published 
luminosity-metallicity (L — Z), mass-metallicity (Af — Z) , 
and mass-effective yield relations. The others provide 
a useful baseline comparison and illustrate some of the 
uncertainties due to choice of calibration and physical 
effects. 



1. iPilvugin fc ThuanI (|2005f) collate recent 
measurements and use them to calibrate 
simultaneous fits to R23 = {/([O II]) -f 
/([O III]AA4959, 5007)}//(H/3) and a func- 
tion of O32 = /([O III]AA4959, 5007)//([O II]). 
The latter is a proxy for ionization parameter, 
which is the ratio of ionizing photons to hydrogen 
nuclei present in gas. 

2. The photoionization models of iMcGaughl ()199lD 
also use both R23 and O32 as parameters in the 
diagnostic. The cali bration is updated and printe d 
in analytic form bv iKuzio de Narav et al.l (|2004f) : 
we use their semi-empirical version. 



3. iTremonti erahl (12004 hereafter T 04) use the mod- 
els of lCharlot fc Longhettil (|200l hereafter GLOl) 
to compute the abundances of '-^10^ galaxies from 
the SDSS. The original models are a suite of ana- 
lytic functions involving different combinations of 
strong emission lines; the choice of diagnostic then 
depe nds on the spectral data available. However, 
ITM directly cross-correlate their data with model 
spectra to find the best-fit abundance. They then 
fit a simple analytic function to the upper branch 
of the abundance vs. R23 relation. 



4. For comparison wit h the | T04| R23 calibration, we 
include the original ICLOlT suite of analytic func- 
tions. We compute abundances using the diagnos- 
tics of Cases A through F, where successive diag- 
nostic cases rely on less spectral information than 
the previous one. For simplicity, we only discuss 
Cases A and F to represent the range of possibil- 
ities in spectral information. The former relies on 
[N II]A6583/[S II]AA6717, 6731 and the latter on 
[O III] A5007/H/3 (with O32 as a weak secondary 
parameter in each case). We note that theiCLQli di- 
agnostics use observed fiuxes as inputs, unlike most 
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Fig. 2. — Abundances in six diflferont diagnostic/calibration systems. (See f\3.1\ for a description of each system.) We select a random 
sample of LIRGs (left) and ULIRGs (right). Colored lines connect different calibrations for the same galaxy. The heavy open circles and 
error bars reflect t he m edian and standard error over the entire sample in each system. We here consider only galaxies that pass the second 
emission-line cut f i|3.3|l . 



other diagnostics which use fluxes corrected for at- 
tenuation by dust. 

5. iKewlev fc Dopital (|2002f ) attempt an optimal cali- 
bration by comparing their photoionization models 
to previous work and combining several diagnos- 
tics. For our galaxies, their 'combined' diagnostic 
reduces to the [N II]A6583/[0 II]AA3727, 3729 di- 
agnostic in almost all cases. 

In each case, we discard galaxies that have log(R23) > 
1; such values are observed in star-forming regions, but 
they are not common (e.g.. iPilvugin fc Thuanll2005r ) . In 
this regime, it is unclear which R23 branch is appropriate 
and the diagnostics are not well-calibrated. Furthermore, 
a very high R23 may be an indication of contributions 
to the line emission from processes not related to star 
formation. There are nine LIRGs and nine ULIRGs with 
R23 > 1. 

We also assume the upper branch for R23-based 
diagnostics ([Edmunds fc Pagei Il984f ). We use the 
[N II] /[O II] flux ratio as a branch indicator, where all 
galaxies with /([N II])//([0 II]) > —1 are assumed to lie 
on the upper branch (e.g., iKuzio de Narav et al.l \20(A 
see also the theo retical plots of [N II]/[ H] and R23 
vs. abundance in lKewlev fc Dopital 12003 ). Eight galax- 
ies (five LIRGs and three ULIRGs) fail this test. Seven 
of these also surpass the R23 threshold; we exclude the 
remaining LIRG from our analysis, to avoid the problem 
of stitching together upper and lower branch calibrations 
(which exist s when two different diagnostics are put to- 
gether; e.g. JSalzer et al...2005) . The number of galaxies 
with low values of [N II] /[O II] is interesting, and may 



suggest an even larger downward spread in abundance 
than is found by taking into account only the upper 
branch systems. 

These two cuts do not strongly impact our analysis. 
Once we put aside galaxies with z > 0.27, we have a 
working sample of 55 LIRGs and 22 ULIRGs. (We dis- 
cuss further the z ~ 0.4 — 0.5 points in §4.41 ) 

Different diagnostic/calibration pairs can yield very 
different abundances for the same galaxy. Figure [2] shows 
the abundance of a random sampling of our galaxies us- 
ing each method listed above. The methods are ordered 
on the horizontal axis by increasing median abundance. 
Over the full sa mple, the extreme media n values (the 
ICLQllCase A and lPilvugin fc ThuanI ()2005f ) cahbrations) 
differ by a factor of 8 for the LIRGs, and correspond to 
12 + log(0/II) = 8.2 and 9.1, respectively. The 'median 
of the medians' is 1.25Zq. For the ULIRGs, the extreme 
values are 8.1 and 9.0, also a spread of a factor of 8, 
with the median being r oughly I.OZq. However, if the 
IPilvugin fc Thuar] ()2005l ) diagnostic (the only one in our 
sample based solely on abundances) is removed, the 
scatter is dramatically reduced. The peak-to-peak varia- 
tion among the five remaining diagnostics is only a factor 
of 2. The median abundances among these five are 1.6Zq 
and IAZq for the LIRGs and ULIRGs, respectively. 

In the next subsections we discuss sources of uncer- 
tainty. 

3.2. Uncertainties in Abundance Caused by Choice of 
Diagnostic 

The lowest abundances arise from the use of the 
empirically-calibrated diagnostics. This exemplifies the 
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12 + log(0/H) [T04] 12 + log(0/H) [C01, case F] 

Fig. 3. — Abundance as a function of calibration for Qeft) [Uharlot & Lonehettil 11200 II) Case F vs. ITremonti et al.l ||2004I '| R23 abundances 
and (right) Chariot & Longhctti (2001) Case A vs. Case F abundances. The former indicates good correspondence on average except at 
high R23 (low iTremonti et a l. 2004 abundances). The latter illustrates good c orrelation but a significa nt offset between the two, independent 
of whether or not the galaxies fall inside the SDSS star-forming galaxy locus IIKauffmann et ani2003l) . Galaxies inside the locus are circled. 
Galaxy designations follow Figure [Tl 



difference between Tg abundances and those from tfieo- 
reticai caiibrations. Di scussi ons of the weaknesses of each 
method are found in iKennicutt et al.l (|2003D and else- 
where, and we do not expand upon them further. How- 
ever, we do note that very recent measurements of tem- 
perature fluctuations in H II regions coupled with heavy 
element recombination line measurements may alleviate 
this d iscrepancy by raising the Te abun dances (iBresolin 
120071) The use of the Te-calibrated iPilvugin fc Thuan 
(|2005( ) method is thus potentially suspect without an ad- 
justment to agree with these findings. Use of their diag- 
nostic calibration is also a particular challenge for LIRGs 
and ULIRGs. We compute low values of their ionization 
parameter variable P ((P) — 0.2) and relatively high val- 
ues of R23 ((log(R23)) = 0.7 for the LIRGs and 0.9 for 
the ULIRGs). These P/R23 combinations are not well 
calibrated for either branch of t he R23 vs. abundance 
diagram in the Pilvugi n fc ThuanI ((20051) diagnostic. 

Of the pho t oioniz ation models, those of 
iKewley fc Dopital (|2002l ) produce the lowest disper- 
sion in our sample, as well as the second-highest 
abundance. The low dispersion is due to the fact that 
the [N II] /[O II] flux ratio is insensitive to the ionization 
parameter. One potential problem with this diagnostic 
calibration is the form assumed for the dependence 
of the N/0 abii n dance ratio on oxygen abundance. 
iKewlev fc Dopital (|2002f ) assume that N/0 is constant 
below 0.23.Z(7i, and linear in abundance above this value 
()van Zee et al.l [T99S' ) . However, the data from recent 
studies seem to show that N/0 may be driven by a linear 
combination of primary and secondary yields, suggesting 
a sum of the two contributions rather than an abrupt 
transition between them. This leads to a slightly higher 
value of N/0 for a given abundance (Bresol in et al.l 
l2005l and references therein). This should lead to a 
higher value of the flux ratio [N II] /[O II] for a given 
abundance, and thus a lower derived abundance from 
the calibration. 



Perhaps the most important calibration issue for the 
result s of this paper is the disagreement between the 
ICLOll Case A— F analytic diagnosti cs and the|TOl empir- 
ical calibration of R23 b ased o n the lCLOll mo dels. Rather 
than comparing all six ICLOll diagnostics to Pfol . we use 
Cases A and F as representative of the greatest and least 
amounts of spectral information required b y the ff LOll 
models. Figures [2] and [3] illustrate that, while 
F and T04 agree reasonably well on average, 
A is higher on average than both by ~0.3 dex. 

What is the cause of this discrepancy? Either the 
LIRGs and ULIRGs have line ratios that are affected by 
processes unrel ated to star formation (see the next sub- 
section) , or the ICLOll Case A diagnostic does not prop- 
erly account for these objects. Looking at Figures [T] and 
H it is tempting t o conclude that Case A, which like the 
Kewlev fc Dopital (j2002) diagnostic is based on a line ra- 
tio involving [N II], predicts higher abundances due to 
[N II] enhancement by a buried AGN or shocks. Con- 
sider the following scenario: AGN and shocks a re imp or- 
tant for some galaxies, for which the different ICLOll di- 
agnostics yield inconsistent abundances; for purely star- 
forming galaxies. Cases A and F instead give the same 
result. More specifically, the objects which fall within 
the SDSS locus of most star-forming galaxies (a crite- 
rion which securely excludes galaxies with even a small 
amount of AGN or shock contribution) sho uld gi ve con- 
sistent results for diagnostics A and F of the lCLOll model, 
while those outside this locus should not. Figure |3| re- 
veals that this is not the case: LIRGs and ULIRGs which 
fall beneath the outer SDSS boundary in the [N II] /Ha 
vs. [O III]/H/3 line ratio diagram (see Fig ure [Tt show 
almost as large a disagreement between the ICLOll Cases 
A and F as do the galaxies that lie outside it. Thus, 
physical effects may have a small contribution, but it ap- 
pears that most of the difference is attributable to the 
inability of the Case A diagnostic to properly account for 
the abundances of these systems. The details of why this 
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Fig. 4.— [N II]A6583/[0 II]AA3727, 3729 vs. 
[O III]A5007/[O II], as measured and predicted by theory. 
These line ratios nicely separate the effects of abundance and 
ionization parameter, and show that starburst models can in 
principle explain LIRG and ULIRG line ratios. Abundance 
increases to the right, with grid points of 0.1, 0.4, 0.8, 2.0, and 4.0 
in units of Zq. Ionization parameter increases upward, with grid 
points of (in the log) roughly 1.2 to 3.0, increasing in steps of 0.3 
dex. The photoionization models are from iKewlev et al. (2001) 
and u se the Starburst99 stellar synthesis code Lcithcrer ct al. 
II1999I ) with 'Lejeune + Schmutz' stellar atmospheres as an input. 
These models assume an electron density of 350 cm~'^ and a 
continuous starburst of age 8 Myr. Galaxy designations follow 
Figure m 

may be so are beyond the scope of this paper. As we dis- 
cuss in the next subsection, other models are consistent 
with all of our LIRGs and ULIRGs being star-forming 
galaxies. 

Among the iCLOll models, the Case F diagnostic mini- 
mizes systematic offsets of LIRGs and ULIRGs wit h re- 
spect to the more sophisticated method by which |T04| 
compute abundances (their method is to find the model 
spectrum that is the best fit to the entire observed spec- 
trum, rather than using just one or two emission-line 
ratios). However, as we see in detail in Figure [3l the 
agreement breaks down at the highest values of R23; in 
this regime, the Case F diagnostic gives higher abun- 
dances by ^0.3 dex. As Figure [2] illustrates, this has 
little effect on the median abundance. Our use oflTOi 
instead of Case F from CLOl thus makes little difference 
to the results of this paper. Despite the fact that it is 
a single-parameter diagnostic, we choose the |T04| diag- 
nostic over the Case F diagnostic because it employs the 
traditional R23 line ratio and is the empirical prod uct of 
a methodologically sophisticated application of the lCLOll 
models. 

3.3. Uncertainties in Abundance Caused by Physical 

Effects 

An obvious important consideration is the dust ex- 
tinction corrections required to derive intrinsic fluxes. 
(The optical attenuation derived from the Balmer decre- 
ment and a foreground screen is F,(B — V) — 0.9 on av- 
erage, corresponding to a factor of 3 upward correction 
to R23 to reach the unattenuated line ratio.) LIRGs and 
ULIRGs have copious quantities of dust in their nuclei; 



the emission from this dust creates their high infrared lu- 
minosities. The optical line emission we observe does not 
emerge from the heavily obscured nuclear star clusters in 
these syste ms, where almost all of the star formation is 
occurrmg (iSurace et al.ll200(D . Presumably it arises in- 
stead from star formation occurring near the surface of 
the nuclear dust cloud, mildly obscured sight lines, or 
light scattered into the line of sight. 

Our use of the Balmer decrement and assumption of 
an effective foreground screen geometry mean that even 
our 'corrected' line fluxes for these dusty systems will 
not correspond to their full line luminosities or star for- 
mation rates. However, several considerations argue that 
the ratios of our corrected line fluxes will still yield fair es- 
timates of abundances. First, ULIRGs in particular show 
such large turbulent velo cities within their compac t nu- 
clear gas structures fe.g.. lDownes fc Solomoiilll998( ) that 
any abundance gradients should be flattened rapidly once 
gas has been dumped into a merger's center of mass. This 
means that the outer and inner regions of a merger's 
dusty core should have similar abundance values. Sec- 
ond, our use of an effective screen geom etry is consis- 
tent w ith both (a) the derivation of the ICalzetti et al.l 
(|2000f ) attenuation law we have adopted, and (b) in many 
objects, the measurement of Balmer decrements much 
higher than could be achi eved if H II regions a nd dust 
were uniformly mixed (e.g.. lThronson et al.lll990h . More- 
over, t ests show that neither us e of a Galactic extinction 
curve (jMiller fc Mathewsl[l97l) nor assumption of a uni- 
formly mixed geometry (when sources' observed Balmer 
decrements allow this) lead to significantly different re- 
sults. Finally, we emphasize that any relative underesti- 
mate of extinction at [O II]AA3727, 3729 would require 
that the true R23 be higher, and therefore the true oxy- 
gen abundance lower, than we have reported here. As 
§§|4] and [5] make clear, such an adjustment would only 
strengthen this paper's conclusions. 

A second relevant issue is the applicability of star- 
formation calibrated abundance diagnostics. As we dis- 
cuss in ij2l LIRGs and ULIRGs deviate from the main 
body of star- forming galaxies in many line-ratio diagrams 
(Fi gure[ll). However, the y generally lie completely below 
the iKewlev et all ()2001f ) maximal starburst lines in at 
least two of these diagrams, and so many or all can plau- 
sibly be treated as starbursts. 

We illustrate this in Figure [4] by plotting the 
[O III]A5007/[O II] line ratio (note this is 0.75 x O32) as 
a function of [N l j]/[0 II] for each gal axy. Photoioniza- 
tion models from iKewlev et ahl (|2001h are plotted atop 
the data, assuming a continuous starburst of age 8 Myr. 
These line ratios nicely separate the effects of abundance 
and ionization parameter. LIRGs and ULIRGs in this 
diagram have relatively high abundances and low ioniza- 
tion parameter, though the details depend somewhat on 
the adopted starburst model. 

It remains possible that some objects in our sam- 
ple may have emission line fluxes with a non-zero con- 
tribution from gas that has been photoionized by an 
AGN or coUisionally ionized by shocks. The addi- 
tion of either a low-luminosity or high-luminosity AGN 
moves galaxies away from the starburs t locu s (e.g., 
IVeilleux fc Osterbrockl ll987: Kcwlcy ct al. 2006b). How- 
ever, it appears that AGN do not make a significant con- 
tribution, because mid-infrared diagnostics of ULIRGs 
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that distinguish starbursts and AGN agree with opti- 
cal diagnostics if all n qn-Seyfert ULIRG s are assumed to 
be star burst-powered (jLutz et al.lll999D . The addition of 
shocks also moves l ine ratios away from the reg ion of pure 
starburst galaxies (Do pita fc Sutherla nd 1995). The ex- 
istence of shocks would be consistent with the presence of 
massive, fast, galaxy-sc ale winds in starburst-dominated 
LIRGs and U LIRGs (iH eckman et al."200lJ: 'Rupke et all 
[2002, 2005bH: iMarthB oOS. 2006). Rupkc ct al. (2005(f 
also find evidence that ULIRGs with LINER spectra 
have higher wind velocities than those with H II spec- 
tral types, suggesting that wind shocks could play a role 
in exciting the observed optical emission lines in LINERs. 

Without careful modeling, we cannot easily flag galax- 
ies with significant contributions from low-level AGN or 
shocks. Rather than attempt an extensive comparison 
to models of stellar and AGN photoionization or shocks, 
we make three empirically and theoretically motivated 
cuts in our data, resulting in three partially overlapping 
subsamples. Computing our results for each subsample 
allows us to ascertain the effect of including galaxies with 
possible low-level AGN or shock contamination. These 
cuts are as follows: (1) accepting on ly galaxies which lie 
beneath the maximal starburst line (jKewlev et al.|[200l[) 
in the [O III]/H/3 vs. [N II] /Ha and [O III]/H/3 vs. 
([S II]AA6717, 6731)/IIa diagrams; (2) accepting only 
galaxies which lie beneath the maximal starburst line 
(|Kewlev et alllMl in the [O III]/H^ vs. [O I]A6300/Ha 
diagram; and (3) accepting only galaxies which lie be- 
neath the line denoting the bound ary of the locu s of 
SDSS star- forming galaxies (Kauffm ann et aLll2003D . As 
shown in later sections, we find that our results are not 
significantly altered by picking one subsample over an- 
other. 

A final uncertainty relates to the apertures covered by 
the spectroscopic observations. Disk galaxies have radial 
abundance gradients, in the sense that abundance de - 
creases with increasing radius fe.g.. lZaritskv et aLlll994D . 
Large apertures covering a significant fraction of the 
galaxy's light thus lead to lower computed abundances 
than those covering only the galaxy's inner regions. For 
the Nearby Field Galaxy Sample (NFGS), covering a 
range of Hubble types, the average difference between 
nuclear and integrated abundanc es is ^0.2 dex for nu- 
clear apertures of 1.5 kpc or less (jKewlev et al.ll200"5l ). 

Our spectra include data with a variety of physical 
apertures, since they are produced with a variety of 
slit widths and extraction apertures. Our LIRG and 
ULIRG samples have median redshifts of 0.04 and 0.14, 
respectively. The LIRG data in particular include one- 
dimensional spectroscopic apertures ~ 1" — 3" in width, 
which corresponds to physical scales of 1 — 2 kpc. 
Th ese qualify more or les s as nuclear spectra according to 
the iKewlev et al.l (|2005D convention for the NFGS. The 
ULIRG apertures also range from ~ 1" — 3" in diame- 
ter, with about half having small apertures. However, 
being at higher redshift, they subtend apertures rang- 
ing from 2 — 7 kpc. Based on comparison to local field 
galaxies, we would then expect a slight systemic devia- 
tion of the average ULIRG to lower abundances relativ e 
to the average LIRG (by ~0.1 dex: lKewlev et al.ll2005f ). 
However, the existence of radial abundance gradients in 
evolved mergers may be inconsistent with expectations of 
strong radial mixing of heavy elements (L. Kewley 2007, 



private communication; lEdmunds fc Greenhowl[r995D . If 
there has been strong radial mixing, aperture size will 
have little affect on the integrated abundance. 

In this work we use two local comparison samples, as 
well as high redshift data. We defer discussions of aper- 
ture effects related to these samples to their points of 
introduction later in the paper. 

With these caveats in mind, we proceed in the next 
sections to employ the internal consisten cy of a given 
abundance diagnostic/calibration pair (theT04R23 cal- 
ibration in our case) to compare LIRGs and ULIRGs to 
other galaxy populations. 

4. NEAR-INFRARED LUMINOSITY-METALLICITY 
RELATION 

The best way to interpret the abundances of local 
LIRGs and ULIRGs in the context of galaxy evolution is 
to compare to the luminosity-metallicity (L — Z) and 
mass-metallicity {M — Z) relationships of other star- 
forming galaxies (i.e., galaxies of lower star formation 
rate or those selected at other wavelengths). These de- 
scribe the increase of abundance with increasing galaxy 
luminosity and mass. They reflect, for a given gas mass 
fraction, either (a) a sequence of changing star forma- 
tion history and enrichment with stellar mass or (b) pref- 
erential accretion of under-enriched gas or loss of over- 
enriched gas in low-mass g alaxies. Recent authors have 
argued for both scenarios jrol : iBrooks et all [20071 : but 
see also iKoppen et al.l l2007l for discussion of an origin 
in a variable initial mass function). We start with the 
L — Z relation, since luminosity is a more easily mea- 
sured quantity than mass. 

We compare to the o nly available near-i nfrared L — Z 
relation for starbursts (jSalzer et al.ll2005[ ). The advan- 
tage of using observations in the near-infrared is their 
lower sensitivity to dust and m ass-to-light ratio vari- 
ations than optic al observations ([Bell fc de Jongj 120001 : 
ISalzer et al.]|2005D . The galaxies comprising this study 
are a subsample of t he KPNO Internat ional Spectro- 
scopic Survey fKISS: [Saker et al.1 [200l . which IS an 
objective prism survey of emission-line gal axies in the 
local universe. The KISS subsample from ISalzer et al.l 
(|2005[ ) consists of galaxies with near-infrared photome- 
try, mostly taken from 2MASS. The galaxies have (z) = 
0.063 and Ks-hand absolute magnitudes ranging from - 
16.5 up to -25 (a factor of 3000 in luminosity). 

Because of the low redshifts and small spectroscopic 
apertures (1.5" — 2", corresponding to 2 kpc) for this 
sample, the computed abundances are approximately nu- 
clear, and we make no aperture corrections. 

4.1. NIR photometry 

Near-infrared photometry for our LIRG and ULIRG 
sample is taken from ground-based data in the i^' band 
where available ([Surace et al.ll2000l : [Stanford et al.|[2000l : 
iKim et all I2002D. and fro m the 2MASS Large Galaxy 
Atlas I Jarrett et al.l |2003[ ) or Extended Source Catalog 
(jCutri et al.ll2006n otherwise. We use 'total,' or extrapo- 
lated, Ks magnitudes from the 2MASS catalog. For some 
objects with multiple nuclei, a 2MASS magnitude is not 
available for the nucleus of in terest but the sky- and star- 
subtracted image is available ijCutri et al.ll2006f ) . In these 
cases, we used the 2MASS Atlas images (in 'postage- 
stamp' form) to do aperture photometry on the nucleus. 
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including associated diffuse emission. In other instances 
a second nucleus was already removed through star sub- 
traction, which was to our advantage. However, in the 
few instances where this subtraction was poor we did 
the photometric separation by hand. For the ULIRGs 
with published nuclear and total magnitudes (Kim e t gdl 
we used the relative nuclear magnitudes to scale 
the total magnitude for the nucleus of interest. Finally, 
we note that the published FIRST-FSC magnitudes are 
on average only 70% of the observed total, as discussed in 
[Stanford et"al1 (|2000.) : in these cases we made an average 
upward correction to reach ^100%. 

We converted measured K' magn itudes to K,^ magni- 
tudes using the formula from Wains coat fc CowTg (jl992D : 
the average H — K' color for H II ULIRGs from 
iSurace et al] (|200l [{H - K ') ^ 0.5, which a pplies 
roughly to LIRGs as well; see lScoville et al.ll2000t ): and 
the color transform ation from CIT sy stem K magnitudes 
to Ks magnitudes (jCutri et al.ll2006l ). 

Starting with the total absolute magnitudes, we con- 
templated a number of adjustments to the measured 
luminosities to reach the rest-frame host galaxy lumi- 
nosity. We decided to apply only one of these ad- 
justments. This involves removal of the central point 
source in ULIRGs, which is most likely an extremely lu- 
minous star c luster (but may a lso contain part of the 
central bulge; ISurace et al.ll2000l ). Very accurate point- 
source subtraction is allowed by high-resolution Hubble 
Space Telescope (HS T) near-infrared data of ULIRGs 
(|Veilleux et al.ll2006l ). The data for five H II galaxies 
from this sample have an average point-spread function 
(PSF) to host-galaxy (PSF-subtracted) luminosity ratio 
of 0.07, with a range of 0.01 to 0.13. Applying this av- 
erage correction to each ULIRG lowers its luminosity by 
an insignificant 0.08 magnitudes. 

We also contemplated corrections for dust extinc- 
tion/emission. Most of the near- infrared light in LIRGs 
and ULIRGs that are dominated by star formation arises 
in modestly extinguished stars, not hot dust. For cold, 
star-forming LIRGs and ULIRGs, most of the near- 
infrared emission i s extra-nuclear (ICarico "eTaLl [l99fll : 
IScoville et all l2000l : IVeilleux et all 12006D, and "the col- 
ors of this extended light are consistent with those of 
normal spirals in the case of LIRGs and a reddened 
stellar population in the case of ULI RGs (C aric o et al.l 
Il990t IScoville et all [2OOOI : iDavies et al. 2002). The av- 



erage predicted screen extinction is modest even for 
the more heavily extinguished nuclear regions (Ay of a 



few; 


Scoville et al.ll2000l: ISurace et al.ll2000l: iDavies et all 


I2OO2 


). To achieve an unreddened luminosity we would 



correct upwar d by only a few tenths of a magnitude 
(jPraind 12003'). and even less for the more dominant 
extended emission. Modest hot dust emission may 
also contribute to the nuclear light (jSurace et al.l 12000 : 
IDavies et al.l l2002l ) , but neither the dust emission nor 
nuclear light is significant. As a final piece of evi- 
dence, ULIRGs tend to show elliptical-like global surface- 
brightness profiles i n the near-infrared, ty pical of relaxed 
stellar populations (jVeilleux et al.l [2OO2I . and references 
therein). 

Finally, we considered K-corrections. There is a 
strong diversity in observed spectral energy distribu- 
tion (SED) shapes in ULIRGs, but the average ULIRG 
shows increasing vLi, from optical to infrared wave- 
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Fig. 5. — Ks band luminosity-metallicity relation for nearby 
emission-line galaxies (small black open circles), LIRGs (blue 
stars), and ULIRGs (red filled circles). Most of the LIRGs and 
ULIRGs fall well below the L — Z relation. The nearby galaxies 
are from the KISS sample, and the black line and dotted li nes are 
a fit t o the data and la RMS dispersion, respectively ( Salz er et al.l 
120051) . The dashed line locates solar abundance. Boxed points do 
not pass our first emission-line cut. We do not plot points with 
z > 0.3. The far-right points do not have measured K magnitudes. 

lengths, and fiat or inc r easing values in the near-infrared 
(jTrentham et al.lll999l : fFarrah et al.ll2003l) . These yield 
positive K-corrections in the if-band that increase slowly 
with increasing redshift and depend somewhat on the 
actual SED shape (0.06 — 0.25 mag at z ^ 0.3 and 
0.1 - 0.4 mag at z - 0.6; iTrentham et al.lll999l ). Given 
the similarity of the near-infrared colors of LIRGs and 
ULIRGs (see previous paragraph and references therein) , 
K-corrections for LIRGs are of similar magnitude. 

We conclude that the K-corrections, extinction correc- 
tions, and corrections for hot dust emission are on aver- 
age small (a few tenths of a magnitude or less in each 
case). Given the inherent variations in SEDs and our in- 
ability to determine precise adjustments for each galaxy, 
we choose to ignore these corrections. Regardless, their 
effect on the following results is negligible. If anything, 
they are likely to increase the significance of the observed 
effects. 

4.2. L-Z relation 

We use the|T03 R23 abundance calibration to compute 
upper branch abundances for our LIRGs and ULIRGs for 
comparison to the KISS results. The resulting A's-band 
L — Z relation is shown in Figures |5| -TTl Each figure 
represents the L — Z relation using a different emission- 
line cut, as described in ^3.3\ It is immediately clear 
that many LIRGs and ULIRGs do not fall on this rela- 
tion, even if we discard possibly shock-excited or AGN- 
contaminated sources. There is some overlap, however. 
We also see that the scatter in the computed LIRG and 
ULIRG abundances for a given magnitude is higher than 
in the KISS comparison sample. 

The observed offset is one of abundance and/or lumi- 
nosity. If the LIRGs and ULIRGs are over-luminous in 
the near-infrared with respect to galaxies of lower star 
formation rates, then there must be a significant (factor 
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Fig. 6. — L - Z relation, but only for LIRGs and ULIRGs that 
pass our second emission-line cut. See Figure [5] for more details. 
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Fig. 7. — L - Z relation, but only for LIRGs and ULIRGs that 
pass our third emission-line cut. See Figurc[5]for more details. 

of ~10) increase in the near-infrared emission above that 
from the old, mass-tracing stellar population. As we dis- 
cuss in t j4.11 hot dust could increase the luminosity above 
that expected from the stellar population, but only by a 
small fraction of the total. Massive red supergiants from 
very young star forming regions can also add to the lu- 
minosity from old stars while keeping the ne ar-infrared 
colors almost constant (|Leitherer et al.|[T999() . However, 
the heaviest young star formation is nuclear, and con- 
tributes only a small percentage of the total near-infrared 
light (SH). 

We conclude that hot dust and young stars cannot 
come close to producing a factor of 10 increase in lumi- 
nosity. Thus, LIRGs and ULIRGs are under-abundant 
with respect to local, emission-line-selected star-forming 
galaxies of lower star formation rate. This conclusion is 
solidified by comparing LIRGs and ULIRGs to the local 
mass-metallicity relation in fj5] 
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Fig. 8. — Difference between the observed abundances in LIRGs 
and ULIRGs and the h — Z relation as a function of infrared lu- 
minosity. The small filled circles represent individual deviations 
from the L — Z relation (from the second emission-line cut). The 
thick open circles are median deviations from L — Z for equal-size 
bins centered on log(Lipi/LQ) = 11.25, 11.75, and 12.25 and the 
error bars represent the standard error in the mean in each bin. 
LIRGs are offset by 0.2 dex, and ULIRGs by 0.4 dex. Comparison 
to the L — Z relation shows a mildly significant trend toward higher 
abundance offsets for higher Lir,. 

The average galaxy in our sample is more NIR- 
luminous than the range of luminosities probed by the 
KISS data (which are biased toward low-luminosity 
objects because o f the required emission-line contrast; 
ISaizer et alll2005l ). Thus, there is uncertainty in the 
shape of the M — Z relation at the luminosities probed 
by our sample. However, even if the relation is conserva- 
tively assumed to flatten above Mk^ ^ —24.5 rather than 
to continue to increase, almost all of the observed galax- 
ies are under-abundant compared to the mean oxygen 
abundance of local, emission-line, modestly star-forming 
galaxies. 

Figure [5] illustrates the offset in abundance from the 
emission-line galaxy relation as a function of total in- 
frared luminosity Lir for the galaxies passing the second 
emission-line cut, and conservatively assuming a flat re- 
lation above Mk^ ~ —24.5. Plotting the data in three 
bins of equal width in log(LiR) shows qualitatively that 
galaxies of higher Lir have higher offsets from L — Z, 
with the LIRGs offset by 0.2 dex and the ULIRGs by 
0.4 dex. Quantitatively, there is a weak and marginally 
significant correlation between the two, with the offset in- 
creasing by ^--^0.2 dex for each 1 dex increase in Lir. The 
correlation coefficient is 0.3 regardless of the emission- 
line cut chosen, while the (parametric) significance of 
the correlation is 99.7%, 98%, and ~90% for the first, 
second, and third cuts, respectively, decreasing largely 
due to the decreasing number of points in each cut. The 
significance of the correlation depends also on the NIR 
luminosity at which the flattening is assumed to occur: 
it increases (decreases) as the cutoff NIR luminosity in- 
creases (decreases). It also depends on the emission- line 
calibration chosen; some calibrations give smaller LIRG- 
to-ULIRG discrepancies (Figure [5]). Finally, we caution 
that the trends may be weaker than suspected from the 
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L — Z relation alone (see Sj5]) . 

4.3. Comparison To Other Merging Systems 

LIRG and ULIRG under-abundances are consis- 
tent with a recent s tudy of optically-selected mergers 
(|Kewlev et al.l [20Q6al) . These authors compute the nu- 
clear abundances of a sample of nearby interacting field 
galaxies with projected separations of 4 — 60 kpc and 
optical luminosities — 22 < < —19. They find an 
under-abundance in mergers with respect to the optical 
i3-band L — Z relation for isolated field galaxies. There 
is a dependence of the offset on separation and 'central 
burst strength,' in the sense that galaxy pairs with small 
separations (4 — 20 kpc) and strong starbursts also have 
the lowest abundances on average compared to isolated 
galaxies. They also find that more optically-luminous 
galaxies have a smaller offset in abundance; for separa- 
tions of 4 — 20 kpc and Mb —21, the abundance offset 
is ^0.1 dex. 

Local LIRGs and ULlRGs have mean optical lumi- 
nositie s comparable t o the highest-luminosity galaxies 

in the iKewlev et all (l2006al) sani ple, with Mr 22 

(|Veilleux et al.l l2002 l: llshid al 120041 ). Based on their nu- 
clear separations (jlshida 20Q2), the non-isolated LIRGs 
may represent similar galaxy interaction states compared 
to the optically-selected galaxies. The mean ratio of 
present-to-past star formation is higher in LIRGs, how- 
ever, based on the higher Ha equivalent w idths in LIRGs 
(|Veilleux et al.l[l995l : iBarton et al.l[200l . The ULIRGs 
represent stronger or more highly progressed galaxy in- 
teractions with much higher star formation rates com- 
pared to the optically-selected mergers. This conclu- 
sion is based on the structural properties of ULIRGs 
(IVeilleux et al.l I2002D and their high Ha equivalent 
widths (jVeilleux et al.lll999| ). More indirectly, ULIRGs 
are scarce in the local universe rela tive to the paren t 
sample of the optically-selected pairs (jFalco et al.lll999l ). 

The observed differences between LIRGs and ULIRGs 
and isolated galaxies of lower star formati on rate reflect 
a continuation of the trends observed by iK ewlev et al.l 
([2006a). LIRGs and ULIRGs are more strongly offset 
from the if -band L — Z relation (by 0.2 — 0.4 dex) than 
are galaxies of similar i?-band luminosity from th e B- 
band L-Z relation (<0.1 dex; IKewlev et al]|2006af) . As 
we describe above, they also are stronger or later stage 
interactions and/or have stronger starbursts. 

Within our sample, the apparent trend of greater abun- 
dance offset with increasing infrared luminosity (Fig- 
ure [8|) mirrors the correlation of abundance offset with 
starburst strength found in optically-selected mergers. 
We also studied the dependence of offset on projected nu- 
clear separation. No significant relationship was found, 
unlike the optical sample. This most likely reflects either 
(a) the saturation of the sensitivity of nuclear separation 
to interaction timescale or merger age as two nuclei merge 
and/or (b) a greater correlation of L — Z offset with level 
of star formation than with nuclear interaction state. 

4.4. Comparison To High Redshift 

Comparison of luminosity-metallicity and mass- 
metallicity relations at low and high redshift point 
to ch emical ev o lution of galaxies over cosmic 
time (iLilly et al. l 120031 iKobulnickv fc KewlevI 
iMaier et all 120051 : ISavaglio et all 120051: lErb et al.l 



2004 



2006 



iMouhcine et al.|[2006D . As mentioned in 21 high-redshift 
abundance measurements of a modest number of LIRGs 
selected at 15 /j,m and ULIRGs selected at submillimeter 
wave l engths exist (Liang et al. 2004; Swinbank et al.l 
12004 iTecza et al.l [2004t iNesvadba et al.l l200iD . Three 
of these s tudies are based on R2.'^ abundances, while 
the third (jSwinbank et al.ll2004f) relies on the [N II] /Ha 
flux ratio, a slightly less robust indicator. From our 
work, we also have data on four LIRGs and ULIRGs at 
z ^ 0.4 — 0.5 (|Rupke et al.ll2005b[ ). Because near-infrared 
photometry for many of these sources exist (iLiang et al.l 
mM iHammer et all [19971: iSmail et all 120041 ). we can 
compare them directly to our data on the if-band L — Z 
relation. 

Of the 15 fim sources from ILiang et al.l ()2004D , we se- 
lect only those which have infrared luminosities consis- 
tent with the definition of LIRGs. We assume the K- 
correction is small and do not app ly it. However, we 
do adjust the magnitudes fisted in (Liang et"al] (|2006l ) 
downwards by 1.9 magnitudes to convert from the AB to 
the Vega pho tometric system , as th ey are listed as AB 
magnitudes in lHammer et al.l (|1997t ). (The magnitude of 
Vega in the AB system is taken from iTokunaga fc: Vaccal 
l2005l ) We also apply upward aperture corrections of 
0.1 dex to reach nuclear abundances, following the treat- 
ment of local disk galaxies in Kewley et al. (2005). 

All of the submi llimeter galaxies (SMGs) in 
ISwinbank et al.l (|2004D have ULIRG-like or higher 
total infrared luminosities. To compare them to the lo- 
cal near-infrared L — Z relation, we apply K-corrcctions 
to the K-hand observed magnitudes from Smail et afl 
(|200l using the avera g e z = 2 ULIRG-derived correction 
from iTrentham et al.l (|1999l ) of 0.45. We also discard 
galaxies with very broad Ha (FWHM > 1500 km s^^) 
and/or log([N II]/Ha) > -0.1. This is an attempt to 
eliminate galaxies with AGN and/or with emission lines 
contaminated by shock excitation, respectively. 

Because of limited emission-line data, we compute 
SMG abundances using the 'coarse' calibration from 
ISalzer et all (|2005f ) that employs [N II]A6583/Ha (and 
is derived from the rT04 R23 relation) . The exception is 
SMM J14011-f0252 for which we use the R23 diagnostic 
with the data from INesvadba et al.l (|2007t l. Comparison 
of the [N II] /Ha and R23 diagnostics for local ULIRGs 
suggest that at low abundances the former tends to over- 
predict the abundance by an average factor of ^--^2. This is 
unsurprising given the results of the application of other 
[N II]-based diagnostics to LIRGs and ULIRGs fi}3?2|l. 

In FigurelH we plot the L — Z data for local LIRGs and 
ULIRGs, for z = OA- 0.9 LIRGs and ULIRGs, and for 
z - 2 SMGs. For SMM J1401H-0252, we plot the abun- 
dances of both nuclei, Jin and Jls. Figure [TUI makes this 
comparison more quantitative. The results show that the 
LIRGs evolve upward in abundance by ^--^0.2 dex from 
z 0.6 to the present day, as would be expected from 
continual processing of heavy elements in their progeni- 
tors. For the z ~ 0.5 ULIRGs, we are limited by number 
statistics, but there is a suggestion of evolution. 

The very high redshift systems, the SMGs aX z ^ 2, 
show a large scatter in abundances. This is partly due to 
low emission-line sensitivity at these redshifts, the coarse 
diagnostic required to compute the abundance, and the 
difficulty in weeding out AGN or shock-excited candi- 
dates. Furthermore, the SMGs are much more near- and 
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Fig. 9. — Redshift evolution of the LIRG (left) and ULIRG (right) L-Z relations. See Figure [5] f or more deta ils. We show here data that 
pass the first emission-line cut. (a) The black, encircled stars are LIRGs with 0.4 < z < 0.9 fromfLTang et aLj (2004, 2006), supplemented 
with one point from this work, (b) The black, encircled circles are ULIRGs with z ~ 0.4 — 0.5 from this work. The black tri angles 
are submillimeter-sel ected galaxies with 1.4 < z < 2.7 and luminosities equal to or greater than those of ULIRGs (Small ct al.i 120041 : 
ISwinbank et al.|[200l) . The vertical line connects the abundances for the two nuclei of SMM J14011+0252 (.Nesvadba et al..,2007.') . 
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Fig. 10. — The abundance off'set from the local L — Z relation for 
low- and high-redshift LIRGs and ULIRGs, as a function of infrared 
luminosity. We show here data that pass the second emission-line 
cut. The (black) small filled circles, (blue) small open diamonds, 
and (red) small open stars represent z ^ 0.1 LIRGs a nd ULIRGs 
from this work, z ~ 0.6 LIRGs and ULIRGs from Liang ct ahj 
(I200I) and this wor k, and 2 ~ 2 SMGs (Swinbank ct al. 2004; 
INesvadba et al.ll2007l ). respectively. The (black) thick open circles, 
(blue) thick open diamonds, and (red) thick open star are median 
deviations from L — Z for local LIRGs and ULIRGs, 0.4 < 2 < 0.9 
LIRGs, and 2 ~ 2 SMGs, respectively. The LIRGs clearly evolve 
upward in abundance by ~0.2 dex from 2 ~ 0.6 to 2 ~ 0.1, as 
would be expected from continual processing of heavy elements. 
Though there only 2 2 ~ 0.5 ULIRGs in this figure, there is also 
apparent redshift evolution in ULIRG abundance from 2 ^ 0.5 to 
2 ~ 0.1. Finally, there is also evidence for modest evolution from 
2 ~ 2 SMGs to 2 ~ 0.1 ULIRGs, though the observed scatter and 
systematic uncertainties are large. 

far-infrared luminous than ULIRGs, so the comparison 
may not be appropriate. Despite these caveats, there is 
a suggestion of abundance evolution on average. 



5. MASS-METALLICITY RELATION 

Galaxy mass is a somewhat more fundamental quantity 
than instantaneous luminosity (i.e., it is more useful for 
predicting a galaxy's properties). The mass-metallicity 
relation is thus a better tool than the luminosity- 
metallicity relation for quantifying the implications of 
LIRG and ULIRG abundances for galaxy evolution. 
However, mass has its own set of limitations; for instance, 
mass is less easily measured than luminosity. In this sec- 
tion, we compare the masses and abundances of LIRGs 
and ULIRGs t o tho se of nearby star-forming galaxies 
from the SDSS (frol l. 

The stellar masses for the ULIRGs are estimated from 
dynamical mass measurements. These are based on near- 
infrared measurements of cent ral stellar velo city disper- 
sions and rotational velocities (jDasvra et al.ll2006bl ) . We 
do not have individual measurements for most (2/3) 
of the ULIRG nuclei in our sample. To plot indi- 
vidual nuclei on the M — Z relation, we instead use 
the high-precision average mass and obse r ved sca tter of 
merger remnant nuclei from lDasvra et all (|2006b[ ) to as- 
sign masses to the ULIRGs in our sample. The variance 
in mass is computed from the mass equation and the ob- 
served scatter in random and rotational velocities. The 
mass of each galaxy is then drawn from a Gaussian ran- 
dom distribution of the proper mean and variance, with 
the limitation that the mass cannot deviate by more than 
3(7 from the mean. 

Direct stellar mass measurements for LIRGs do not 
exist in the literature. However, the indirect method 
of modeling galaxy spectra suggests an average stellar 
mass in star-forming LIRGs that is one - half that of a 
ULIR G merger remnant (^ 5 x 10^° M^ HPasquah et all 
I2OO5I ). This is co nsistent with some LIRGs being pro- 
genitor ULIRGs (llshidal |2004 ). A caveat is that the 
iPasauah et al.l ()2005f ) galaxies have log(LFiR/i0) = 
10.5— 11.5, rather than the range log{Lm/ Lq) — 11.0 — 
12.0 normally assigned to LIRGs (note that Lpm < Ljr 
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Fig. 11. — Comparison of the mass-metallicity relation from the 
SDSS l ITremonti et al.ll200l) with LIRG and ULIRG abundances 
and stellar masses. The average LIRG and ULIRG are significantly 
under-abundant, as they are when compared to the L — Z relation. 
The dotted lines show Icr scatter on either side of the mean SDSS 
relation, which has been shifted upward by 0.1 dex to account for 
aperture effects. Atop this relation are median LIRG and ULIRG 
abundances (colored diamonds) . The diamond colors represent dif- 
ferent emission-line cuts (cut 1, green; cut 2, blue; cut 3, black), 
and the sizes represent the dispersion in points (the standard error 
in the median). We also plot individual abundance measurements 
under the second emission-line cut, randomly distributed in mass 
according to t he measured m ean and standard deviation for LIRGs 
and ULIRGs llPasguali et al. 2005; Dasyra ct al. 2006b). 



by ~10%; ISanders fc Mirabell [TOOl) . Regardless, our 
comparison to M — Z does not depend sensitively on our 
choice of LIRG mass. As with the ULIRGs, the mass of 
each galaxy is drawn from a Gaussian random distribu- 
tion of the proper mean and variance. We assume that 
the mean equals one-half the mean for ULIRGs but that 
the variance in linear mass space is the same. 

The SDSS galaxies to which we compare are primarily 
late- type (TO^), so their disks presumably have abun- 
dance gradients. The physical diameter spanned by the 
SDSS spectroscopic fiber is relatively large (3", corre- 
sponding to 4.6 kpc at (z) — 0.08). Using trends of 
abundance dilution as a function of aperture size gleaned 
from the NFGS data (Kcwley et al. 2005), we correct the 
SDSS abundances upward by 0.1 dex so that they better 
approximate nuclear abundances. 

The M — Z relation is shown in Figure [TT] with the 
LIRG and ULIRG points over-plotted. The results are 
comparable to those obtained when we compare LIRGs 
and ULIRGs to the L - Z relation. LIRGs and ULIRGs 
are significantly offset from the M—Z relation, regardless 
of the emission-line cut chosen. We also observe a much 
larger scatter in the abundances of LIRGs and ULIRGs 
than in the reference sample at a similar mass. 

In the case of M — Z, the offsets from the mean rela- 
tion are unambiguously ones of abundance rather than 
mass, since the progenitors of LIRGs and ULIRGs are 
massive spirals. In Figure [T^ we plot the abundance off- 
sets of LIRGs and ULIRGs from the M — Z relation as 
a function of infrared luminosity atop the same offsets 
computed using the L — Z relation. In each case we use 
data that pass the second emission-line cut. The median 
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Fig. 12. — Difference between the observed abundances in LIRGs 
and ULIRGs and the L — Z and M — Z relations, as a function 
of infrared luminosity. The (red) thick open diamonds are median 
deviations from the M — Z relation. There is excellent qualitative 
agreement between the two, though comparison to the M — Z re- 
lation yields a slightly larger under-abundance of the LIRGs than 
does comparison to the L — Z relation. Refer to Figure [8] for more 
details. 

offset of ULIRGs from the M — Z relation, 0.4 dex, is 
the same as the offset from the L — Z relation (with our 
conservative flattening assumption), but the offset from 
the M — Z relation is higher for the LIRGs at 0.3 dex 
(vs. 0.2 dex from the L — Z relation). The correlation 
between Xir and under-abundance may be present in 
the M — Z offsets as it is in the L — Z offsets, but the 
observed trend is weaker and not significant. Neverthe- 
less, given the systematic uncertainties in abundance di- 
agnostic calibrations, individual luminosity corrections, 
individual masses, and the exact shape of the L — Z re- 
lation, we conclude that the L — Z and M — Z offsets are 
in good agreement with one another. 

6. EFFECTIVE YIELD 

Using our measurements of the gas-phase abundances 
of oxygen in LIRGs and ULIRGs, we are able to com- 
pute effective yields. The true yield p is the fraction 
of the mass of a generation of stars that is converted 
into a heavy element (in this case, oxygen) and returned 
to the ISM. More precisely, for a given stellar gener- 
ation, -p refers to the total mass of a heavy element 
produced by massive, short-lived stars normalized by 
the mass locked up in long-lived stars and stellar rem- 
nants. The related quantity of effective yield is defined 
as pe// = Z I ln{\i-^), where = Mgo.sl{Mgas+MstaTs\ 
is the gas mass fraction and Z = Mheavy element /M gas- 
The effective yield provides information on the chemical 
history of the galaxy through comparison with detailed 
evolutionary models. In the case of a 'closed-box' model 
with instantaneous recycling, the effective yield equals 
the true yield {p = Peff)- 

The effective yield is more sensitive to the chemical his- 
tory of galaxies than the M — Z relation alone, since it 
also incorporates information about the present gas con- 
tent of the galaxy. Star formation increases a galaxy's ef- 
fective yield until it asymptotically reaches the true yield. 
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by consuming gas and producing metals. Conversely, 
gas flows in and out of the galaxy reduce the ef fective 
yield (e.g lEdmundsl Il990t iKoppen fc Edmunds! Il999t 
iDalcantonI |2007[ ). An exception is the case of inflows of 
gas with non-zero abundance. Elementary models of ra- 
dial inflow in galaxy disks suggest that they can increase 
the effective yield a bove the true yield by small a mounts 
(factors of 2 or less; lEdmunds fc Greenhow|[l995l ). 

A typical ULIRG gas mass fraction can be deter- 
mined from global H I and H2 observations. Some un- 
certainty exists due to the fact that the conversion of 
CO molecular luminosities to H2 gas masses differs be- 
tween UL IRGs and norm al gala xies and is not known ex- 
actly fcf. JSanders et al.' 1991; D ownes &: Solom on 199^. 
The uncertainty is compounded by the scarcity of H I 
measurements of ULIRGs. We use the best available 
data on ULIRG gas masses (iMirabel fc SandersI fl988L 
ll989HSanders et al.lll99llDownes fc Solomonlll998D . and 
the average s t ellar m ass of single nucleus ULIRGs from 
iDasvra et al.l (|2006af l. We then consider ranges of possi- 
ble values for the H I mass fraction and CO-to-Il2 con- 
version factor, arriving at an average gas mass fraction 
in ULIRGs of (/ig) = 0.1 ± 0.05. (The conservative un- 
certainty estimate is dominated by our uncertainty in 
Mhi IMh2 ■ ) This estimate is roughly consistent with the 
gas mass fraction estimated fr om gas dynamical measure - 
ments (Mgas/Mdyn ^ 0.16; iDownes fc SolomonI fl998h . 
As expected, it is also much lower than the gas mass frac- 
tions of high-re dshift, subniillimeter-selected ULIRG s 
iiHg) = 0.3-0.5: [C^ve et al.l[2005t iTacconi et al.ll2006n . 

In the case of LIRGs, Mhi/Mh-, is bette r con- 
strained by measu rements ()Mirabel fc Sanderi 119891 : 
ISanders et al.l fl99lb . The value of Mj/^Z-^co ^^^o 
lower than in normal galaxies, though its value relative 
to that in ULIRGs is uncertain. An additional uncer- 
tainty arises due to the absence of accurate stellar mass 
measurements in LIRGs. We resort to assuming that 
the stellar mass of a LIRG is one-half that of a ULIRG 
(as in !j5]). Again, using the best measurements avail- 
able and bracketing different possibilities, we estimate 
(Hg) = 0.2 ± 0.1. (In this case, the conservative uncer- 
tainty estimate is dominated by our uncertainties in the 
stellar mass and MH2/ L'qq.) 

Figure [13] displays the median effective yie lds o f oxy- 
gen m LIRGs and ULIRGs compared to the ffp mass- 
effective yield relation. LIRGs and ULIRGs have sig- 
nificantly lower effective yields, by factors of 2 and 3.5, 
respectively. This result mirrors the under-abundances 
described in previous sections, since (a) Peff is linearly 
proportional to oxygen mass fraction and (b) the LIRG 
and ULIRG gas mass fractions (0.2 and 0.1, respectively) 
are similar to those of the average SDSS emission-line 
galaxy of similar mass (/Xg = 0.2 for M ~ 10^^ Mq). 
This result is robust even if we have seriously underes- 
timated jjLg in ULIRGs: if ^g in ULIRGs is raised to 
0.2, pe// is still significantly lower (by a factor of 2.5) in 
ULIRGs than in the SDSS galaxies. 

7. DISCUSSION 

What is the cause of the L — Z and M — Z under- 
abundances (Figures[5]-l7landfTT|l? It is accepted that the 
average progenitors of LIRGs and ULIRGs are gas-rich 
spira l s (iBarnes fc Hernquis^ 119961: iMihos fc Hernguistj 
I199I I Veilleux et al.ll200 J 1 lono et al. 1 120041 : llshidal 120041 : 
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Fig. 13. — LIRG and ULIRG average effective yields compared 
to th e mass-effective yield relation from the SDSS l ITremonti et al.l 
120041) . The average LIRG and ULIRG have effective yields that are 
lower than the average SDSS galaxy of the same masses by factors 
of 2 and 3, respectively, as a consequence of abundance dilution. 
The yields are computed using the median abundances for LIRGs 
and ULIRGs and gas mass fraction esti mat es of fig = 0.1 ± 0.05 
and 0.2 it 0.1, respectively. See Figure [TT] for an explanation of 
most of the lines and symbols. The horizontal dashed line is the 
SDSS estimate of the true yield. The SDSS yields have been shifted 
upward by 0.1 dex to account for aperture effects. 



INaab et al.l[200l . Prior to an interaction and the en- 
suing intense star formation episode, these progenitor 
spirals presumably lie on the L — Z and M — Z re- 
lations. Since star formation increases the abundance 
of a system, the only proce sses that can reduce it are 
gas motions (lEdmundslll990HKoppen fc Edmunds! 1 1999t 
!Dalcantonl!2007[) . 

Almost all LIRGs and UL IRGs ho s t fast , massive, and 
powe r ful ou tfl ows ( Hcckm an et al.! !2000! : ! Rupke et al.l 
[2OOI [2005aH : IMartinl 120051 . !2006!) . The hot phase of 
these outflows is en riched in oxygen relative to iron (e.g., 
!Grimes et al.|[2005l) . This enrichment arises because the 
hot phase consists of supernovae ejecta from the current 
generation of star formation. Expulsion of these metal- 
rich ejecta therefore cannot in principle reduce the abun- 
dance of the galaxy below its value prior to the starburst. 
Rather, outflow of hot gas prevents further enrichment. 

There exists also a neutral, much cooler gas phase in 
these o utflows with mass s everal times that in the hot 
phase (jRupke et al.l !20053 )^. The average neutral gas 
mass in LIRG and ULIRG winds is estimated to be 
10^ Mq (IRupke et al l 120053) . or ~10% of the average 
total gas mass. The effect on the galaxy's ISM is thus 
significant, and outflows could remove the bulk of the 
galaxy's ISM over long timescales (10^-10^ yr) if they 
act continuously. The outflowing neutral gas is swept 
up from the disk of the galaxy and has abundance equal 
to that of the nuclear star-forming region. Therefore, 



^ IRupke et al.l Il2005d) write that the mean M^aldl ^'^hot in 
LIRGs and ULIRGs is ~0.5. However, the normalization of the 
model that was compared to the observations is in error. Correct- 
ing it implies that the actual entrainment value is ~10 in LIRGs 
and ~2 in ULIRGs, with an average of a few and a large scatter 
(from 0.25 to 20). 
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despite its large mass, the neutral phase of LIRG and 
ULIRG outflows cannot reduce the abundance of the nu- 
clear regions, since it is not preferentially removing heavy 
metals from the ambient gas. 

Significant radial inflo w of gas also occ urs in 

gas- rich mergers (e.g. , iBarnes fc HernauistI 119961: 
Mihos fc HernauistI 119961: llono et all 12004 iNaab et all 



20061) . Transfer of angular momentum from gas to 



stars through gravitational torques and dissipation of 
kinetic energy in shocks causes much of the gas of the 
progenitor s ystem s to flow inward during the merger. 
iNaai] et "aTl (|2006( l study the merger of two disk galaxies 
of equal mass, each with /ig = 0.1. On the galaxies' 
first pass (the LIRG phase), 10 — 60% of the gas in each 
disk flows inward to small radii (with the exact value 
depending on geometry). During the final coalescence 
(the ULIRG phase), 50 — 85% of the total gas moves to 
the nucleus. The final increase in the centr al gas mass 
is by a factor 3 — 4 within a 2 kpc radius (jlono et al.l 
[200l . 

Radial abundance gradients exist in most disk galaxies, 
such that abundances of oxygen and other elements de- 
crease MTth_^ncreasmg_dls from the galactic center 
(e.g.. IZaritskv et aLlll994D . Inflow of gas consequently 
moves gas of low abundance into a nuclear region whose 
initial abundance is higher. As a result, the central 
regions of LIRGs and ULIRGs would naturally exhibit 
lowered abundances due to dilution. The dilution of 
the nuclear abundance Z is given by Z final I Z initial = 
(1 -f a6)/(l -I- a), where a is the ratio of gas mass ar- 
riving in the nucleus to that of gas already present and 
h is the ratio of abundances of the inflowing and nu- 
clear gas. For LIRGs and ULIRGs we have shown that 
Z final I Z initial ^ 0.5, and simulations predict a = 3 — 4 
(|Iono et al.l2004f) . Inverting the above equation and solv- 
ing for 6, we find h = 0.35 (i.e., the inflowing gas had 
35% of the initial abundance of the nuclear gas). This 
is consistent with the mean oxygen abundance gradi- 
ent observed in spirals (-0.6 dex per isophotal radius, 
or 25% of the nuclea r abundance at the isophotal radius; 
IZaritskv et"al]|1994f ). 

The effective yield of a galaxy can also change under 
a variety of processes (§[!]). However, in our case the 
low effective yields we observe are primarily due to low 
abundances. There is a secondary effect on Peff due to 
changes in /ig, but the effect is small for these systems 
(®. Thus, the effective yields in LIRGs and ULIRGs 
have been lowered primarily by the same processes that 
have lowered the abundances. 

The inflow scenario was proposed to explain the off- 
set of local, optically-sel ected interacting pair s from the 
i?-band L — Z relation (|Kewlev et al.l [2C)06a') . Our ar- 
guments above, based on reduced abundance and effec- 
tive yield, clearly make it the preferred explanation for 
the se in frared-selected mergers, as well. As we discuss 
in ij4.3[ the offsets from the L — Z relation in LIRGs 
and ULIRGs are larger than those in optically-selected 
mergers. Our interpretation implies that more gas has 
been funneled to the centers of LIRGs and ULIRGs than 
in the optical mergers. This parallels the observation 
that ULIRGs are later-stage mergers and that LIRGs 
and ULIRGs have higher star formation rates than the 
optical mergers. The amount of gas channeled to th e 
center increases with merger age (e.g.. lNaab et al.ll2006f ). 



and star formation is powered by the gas dr iven to the 
merger center (e.g.. lMihos fc Hernquistlll996l ). 

The differences between the optical and infrared merg- 
ers could also arise from their different chemical evolu- 
tionary histories following inflow. We have concluded 
that under-abundances in the infrared-selected systems 
must be caused by gas inflow. Following a major in- 
flow event, closed-box consumption of gas by star for- 
mation in LIRGs and ULIRGs will tend to elevate the 
gas-phase abundances of heavy elements, while strong 
outflows of metal-rich supernovae ejecta will maintain 
low abundances. Thus, the abundance dilution due to 
inflow may have initially been stronger than observed, 
depending on the subsequent effects of star formation 
and outflow. 

The degree to which the initial under-abundance cre- 
ated by inflow was smaller than the observed value can 
in principle be determined from the properties of the en- 
suing star formation and outflow. In practice, the degen- 
eracy among quantities such as the amount of gas con- 
sumed by star formation, the amount of gas ejected by 
outflows, and the enrichment of the outflowing gas means 
that such an exercise is not well-constrained. Simply 
from the good agreement between the observed under- 
abundances and those expected from our basic calcula- 
tions of abundance dilution, we tentatively conclude that 
the abundance immediately after inflow cannot be much 
lower than the abundance observed now. In Appendix El 
we present example calculations that show subsequent 
abundance evolution due to star formation and outflows 
need not undermine this argument. The uncertainties 
are too large to make quantitatively useful statements 
about the impact of star formation and outflows. How- 
ever, these calculations do serve to illustrate that the 
scatter in abundances we observe may be due to galaxy- 
to-galaxy differences in how the inflow, star formation, 
and outflow have proceeded. 

8. SUMMARY AND OUTLOOK 

We have shown that luminous and ultraluminous in- 
frared galaxies do not follow the standard luminosity- 
metallicity and mass-metallicity relations for isolated 
galaxies with lower star formation rates. As a conse- 
quence, their effective yields are also smaller than those 
of similar-mass galaxies. 

We conclude that there is an oxygen under-abundance 
in the nuclei of LIRGs and ULIRGs, compared to galaxies 
of the same luminosity and mass. We attribute this to ra- 
dial inflow of gas into the galaxy nuclei, as was concluded 
for studies of optically-selected mergers (jKewlev et al.l 
[2006a). This is consistent with the fact that many LIRGs 
and all ULIRGs are in the early and late stages, respec- 
tively, of a merger of two roughly equal-mass galaxies. 
Other LIRGs are involved in more minor interactions 
and/or mergers. These interactions and mergers cause 
gas from large radii, which is less oxygen-abundant than 
the central regions, to fall to the center and dilute the 
nuclear abundance. This infalling gas also fuels star for- 
mation; LIRGs and ULIRGs manifest the highest star 
formation rates in the local universe. 

Consistent with this interpretation is the observation 
that the LIRGs and ULIRGs have larger offsets from the 
L—Z relation (0.2—0.4 dex) than o ptically-selected merg - 
ers of similar luminosity (<0.1 dex: lKewlev et aLll2006af) . 
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The LIRGs and ULIRGs are in more extreme interaction 
states and/or have much higher star formation rates than 
the optically-selected mergers. 

A consequence of this result is that the L — Z and 
M — Z relations are not universal. At z > 0.5, where 
the star formation rate density of the universe is domi- 
nated by LIRGs an d ULIRGs fe.g.- lLe Floc'h et al.ll2005l : 
ICaputi et all 120071) . the L - Z and M - Z relations 
for star-forrning g a laxies wil l be governed by rnergers 
(jSmail et all |2004; Shi et al.1 [200a iBridee et al.l l2007l . 
though see also Melb ourne et al.ll2005( ). Observations of 
high-redshift L — Z and M — Z relations that reveal a 
relation shifted to lower abundance and/or higher lumi- 
nosity and mass must account for this. For instance, the 
optically-selected galaxies at z 2 that were use d to as- 
sess the evolution of M — Z are mostly LIRGs (,Erb et al.l 
120061: iReddv et al1l2OO60 . A large part of the offset at- 
tributed to evolution may be due to the fact that the re- 
lations include more galaxies that are the by-product of 
merging. These high-z mergers would have lower abun- 
dances not only because they are younger but also be- 
cause they have undergone gas inflow. 

Our low-redshift sample thus represents a baseline for 
comparison to abundances of high-redshift LIRGs and 
ULIRGs. The abundances of high-redshift luminous in- 
frared galaxies should be compared only with the abun- 
dances of local galaxies with comparable star formation 
rate and/or interaction strength, not with the abun- 
dances of field galaxies of low star formation rate. We 
compare our sam ple to a small sample of LIRGs at 
0.4 < z < 0.9 (Liang et al.l I2004D and find an increase 
in abundance with decreasing redshift (~0.2 dex from 
z - 0.6 to z = 0.1). Three z ~ 0.5 ULIRGs from 
the current sample also have lower abundances than the 
z = 0.1 mean. The situat ion for ex tremely lu minous, 
high-redshift SMGs^ (Swin bank et al, 2004; Tccz a et al.1 
l2004t iNesvadba et alliSlOOTf ) is unclear, however; our up- 
per limit on the mean abundance is consistent with red- 



shift evolution but more data are needed. Many more 
optical and near-infrared spectra of z > 0.5 LIRGs and 
ULIRGs are necessary to understand the abundance evo- 
lution of infrared-selected mergers. 
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APPENDIX 

MODELS OF POST-INFLOW ENRICHMENT 

To illustrate interstellar chemical evolution in LIRGs and ULIRGs subsequent to gas inflow, we present two example 
calculations, one using a simple closed-box and the other a modified leaky-box model. In these examples, we assume 
the true oxygen yield of p = 0.0131 from iT04 . adjusted upward by 0.1 dex according to our aperture corrections. We 
refer to the oxygen abundance by mass after inflow and before star formation and outfl ow as Zi, and the observed 
abundance as Z2- We use the average observed value of Z2 ~ 0.006 from this paper (in the |T04| calibration). Similarly, 
the galaxy's gas masses before and after star formation/outflow are Mi and M2. For simplicity, we choose a single 
value of M2/M1 = 0.8 as a reasonable guess (i.e., 20% of the galaxy's gas has been consumed in star formation and/or 
ejected by outflows). Exactly how much gas has been consum ed by star formation is u nclear, but we do know that 
the gas-consumption timescale in ULIRGs is short (<100 Mvr: lDownes fc Solomonlll998l ). Because our goal is simply 
to illustrate that Z2/Z1 ca n be minimized , we leave it to5B the reader to try other values of M2/M1. 

In the closed-box model (jTalbot &: Arn ett 1971), Zi is given by 

Zi (closed box) = Z2+P ln(M2/Afi). (Al) 

Substituting the values listed above, we find Zi = 0.003, or half the present-day value. In the closed-box model with 
the values listed above, consuming more than half of the gas present (M2/M1 < 0.5) implies that Zi < 0, which is 

unphysical. 

In a leaky-box model (jHartwickl 119761 ) . outflows eject from the galaxy some or all of the heavy elements produced 
by ongoing star formation. Two new quantities enter: rj = ('^^J^^^of mass outflow rate normalized by the star 

formation rate, and a = Zpp / Z jsm is the enrich ment of the outflow relative to the ambient ISM . For our derivation, 
we begin with the treatment of iMatteuccil (|200lD . We modify equation 5.33 of iMatteuccil (|2001f ) by adding a to the 
third term: 

d{ZMg) /dt ^ p SFR - Z SFR -a Z {dM/dt)oF, (A2) 
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where Z is the ISM abundance, Mg is the gas mass, SFR is the star formation rate, and [dM j d£)oF is the mass 
outflow rate. All quantities but p, 77, and a are assum ed to be a function of time. Using the fact that dMg/dt — 
— {dM/dt)oF — SFR (equation 5.32 of lMatteuccill200ll ). we solve this equation through substitution and integration. 
The result is the equation for the modified leaky-box case of an outflow with abundance greater than that of the 
ambient ISM: 

Zi (leaky box, a>l)=p/c- {Ah/Ahy'/'^ip/c - Z2), (A3) 

where c = r]{a — 1) and d = i] + 1. Because star formation and outflows cannot lower the galaxy's abundance (i.e., 
Zi < Z2; see 33)7 we have the further constraint that p/c — Z2 > 0, or c < P/Z2 (c < 2.18, for our values of p and Z2). 
In the case a = 1, the proper equation is the standard leaky-box case of an outflow with abundance equal to that of 
the ambient ISM: 

Zi (leaky box, a ^ 1) ^ Z2 + {p/d) ln(M2/A/i). (A4) 

iRupke et al.l (|2005cD estimate {rf) ~ 0.2 — 0.3 (with a large scatter) for the cool outflow phase of LIRGs and ULIRGs. 
The cool phase, which is mostly entrained ambient ISM, has Zqf ~ Zjsm, or a = 1. For star formation plus a 
cool-gas outflow, Zi = 0.0036 — 0.0039, which is not far from the closed-box case. Determining 77 and a for the hot 
phase of LIRG and ULIRG outflows is quite difficult, as they are poorly con strained obse r yation ally. By comparing 
the mass outflow rate in the cold phase to predictions from starburst models, iRupke et al.l (|2005d ) find 77 could be as 
high as ^0.1. Through m odel fits to the a-element enhancement. X-ray observations suggest a could reach values of 
-10 dG rimes et al.ll2005l ). These values yield Zi = 0.0045, 25% lower than the observed abundance. If we combine 
the hot and cold phases, we have ij ~ 0.3 — 0.4 and (a) ~ 4 — 4.3, which yields Zi = 0.0046 — 0.0051. This is only 
15 — 20% lower than the average observed abundance in LIRGs and ULIRGs. 

This exercise, though difficult because of observational uncertainties, demonstrates two things: (a) the amount of 
post-infiow enrichment required to reach the observed values of oxygen abundance in LIRGs and ULIRGs need not be 
large, and (b) galaxy-to-galaxy variations in chemical evolutionary histories following infiow events can easily produce 
part of the observed scatter in abundances (through, e.g., varying degrees of gas consumption, outflow enrichment, 
and outflow efficiency). 
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TABLE 1 
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F00267+3016:NW 


2MASXi J0029256+303325 


0, 


,050 


11 


.73 


-25, 


,39 


76 


.80 


2jy 


3 


F00267+3016:SE 


2MASX J00292498+3033339 


0, 


,050 


11 


.21 


-24, 


,09 


23 


.19 


2jy 


3 


F01173+1405 


2MASX J01200265+1421417 


0, 


,0312 


11 


.63 


-24, 


,38 


100 


.00 


rbgs 


423 


F01364— 1042 


2MASX J01385289— 1027113 


0, 


,0484 


11 


.76 


-23, 


,88 


100, 


.00 


rbgs 


4 


F01484+2220 


2MASX J01511437+2234561 


0, 


,0323 


11 


.64 


-25, 


,50 


100, 


.00 


rbgs 


5 


F02248+2621 


2MASX J02274641+2635222 


0, 


,033 


11 


.42 


-24, 


,59 


100 


.00 


2jy 


3 


F02512+1446 


UGC 2369 NEDOl^ 


0, 


,0315 


11 


.70 


-24, 


,36 


100, 


.00 


rbgs 


31 


Fl_5 


[CDF99] Fl_005 


0, 


,4786 


11 


.85 


-1 


,00 


100, 


.00 


cdf 


1 


F04315— 0840 


2MASX J04340002— 0834445 


0, 


,0160 


11 


.60 


-24, 


,74 


100, 


.00 


rbgs 


6 


F06538+4628:SW 


2MASX J06573445+4624108 


0, 


,021 


11 


.31 


-24, 


,03 


100, 


.00 


rbgs 


3 


IRAS07062+2041 


2MASX J07091189+2036102 


0, 


,017 


11, 


.13 


-24, 


,17 


100, 


.00 


2jy 


3 


IRAS07063+2043 


2MASX J07091808+2038092 


0, 


,017 


11, 


.26 


-24, 


,80 


100, 


.00 


rbgs 


3 


F07256+3355 


2MASX J07285341+3349084 


0, 


,013 


11, 


.23 


-24, 


,16 


100 


.00 


rbgs 


3 


F08572+3915:SE 


IRAS 08572+3915NW 


0, 


.0583 


11, 


.23 


-22. 


.17 


13, 


.57 


ijy 


4 


F09046+1838 


2MASX J09073082+1826057 


0. 


,029 


11, 


.31 


-24. 


.78 


100, 


.00 


2jy 


3 


F09126+4432:NE 


2MASX J09155548+4419576 


0, 


.0398 


11, 


.70 


-25. 


.28 


100, 


.00 


rbgs 


4 


F09209+3943:E 


2MASX J09240034+3930426 


0, 


,0921 


11, 


.52 


-25, 


.17 


65 


.29 


wgs 


24 


F09209+3943:W 


2MASX ,109235974+3930456 


0, 


,0922 


11 


.25 


-24, 


,49 


34, 


.70 


wgs 


4 


F09218+3428 


FIRST J092455. 0+341535 


0, 


,068 


11 


.72 


-1 


,00 


100, 


.00 


wgs 


2 


F09320+6134 


2MASX ,109355169+6121105 


0, 


,0393 


11, 


.96 


-25, 


,53 


100, 


.00 


rbgs 


423 


F09333+4841 


2MASX J09363719+4828275 


0, 


,0259 


11, 


.37 


-24, 


,20 


100, 


.00 


rbgs 


3 


F09338+3133 


2MASX J09364832+3119507 


0, 


,077 


11 


.62 


-24, 


,86 


100 


.00 


wgs 


2 


F09339+2835 


2MASX J09365478+2822200 


0, 


,119 


11 


.78 


-24, 


,96 


100, 


.00 


wgs 


2 


F09399+2830 


2MASX J09425370+2816563 


0, 


,053 


11 


.17 


-23, 


,58 


100, 


.00 


wgs 


2 


F10190+1322:W 


F10190+1322W 


0, 


,0766 


11 


.68 


-24, 


,93 


41 


.56 


Ijy 


1 


F10190+1322:E 


F10190+1322E 


0, 


,0759 


11 


.82 


-25, 


,27 


58, 


.44 


Ijy 


41 


F10203+5235 


2MASX J10233258+5220308 


0, 


,0322 


11 


.43 


-24, 


,94 


100, 


.00 


2jy 


34 


F10594+3818:SW 


F10594+3818SW 


0, 


,1578 


11 


.94 


-24, 


,68 


38, 


.03 


Ijy 


4 


F11231+1456 


2MASX ,111254505+1410359 


0, 


,034 


11 


.57 


-24, 


,81 


100, 


.00 


rbgs 


3 


F12116+5448 


2MASX J12140957+5431360 


0, 


,0081 


11, 


.06 


-23, 


,09 


100, 


.00 


rbgs 


64 


F12251+4026 


2MASX J12273800+4009378 


0. 


,037 


11, 


.49 


-24. 


.72 


100, 


.00 


2jy 


3 


F12592+0436 


2MASX J13015026+0420005 


0. 


,0373 


11, 


.62 


-24. 


.84 


100, 


.00 


rbgs 


4 


F13136+6223 


2MASX J13153506+6207287 


0. 


,031 


11, 


.79 


-24. 


.18 


100, 


.00 


rbgs 


4236 


F13182+3424 


2MASX J13203537+3408218 


0. 


,022 


11, 


.67 


-23. 


.98 


100, 


.00 


rbgs 


56 


F13373+0105:NW 


2MASX J13395227+0050224 


0. 


.022 


11, 


.39 


-24. 


.82 


100, 


.00 


rbgs 


6 


F13373+0105:SE 


2MASX J13395767+0049514 


0. 


,023 


11, 


.15 


-25. 


.04 


100 


.00 


rbgs 


6 


F13458+1540:NE 


2MASX J13481477+1525497 


0, 


,058 


11 


.43 


-24, 


,96 


49, 


.07 


2jy 


3 


F13458+1540:SW 


CGCG 102—061 NEDOl^ 


0, 


,058 


11 


.45 


-25, 


,00 


50, 


.92 


2jy 


3 


F14179+4927 


2MASX J14194323+4914121 


0, 


,026 


11 


.32 


-23, 


,97 


100, 


.00 


rbgs 


5 


F14280+3126 


2MASX J14301041+3112558 


0, 


,012 


11 


.06 


-24, 


,40 


100, 


.00 


rbgs 


5 


F14416+6618 


2MASX J14423488+6606043 


0, 


,038 


11 


.27 


-23, 


,38 


100, 


.00 


wgs 


2 


F15043+5754:S 


F15043+5754S 


0, 


,1505 


11 


.94 


-25, 


,44 


49, 


.54 


Ijy 


4 


F15163+4255:NW 


2MASX J15180612+4244445 


0, 


,039 


11 


.90 


-24, 


,63 


100, 


.00 


rbgs 


5 


F15359+3139 


2MASX J15375657+3129582 


0, 


,0534 


11 


.02 


-23, 


,82 


100, 


.00 


wgs 


24 


F15384+3841 


2MASX J15401502+3832111 


0, 


,0673 


11 


.24 


-24, 


,82 


100, 


.00 


wgs 


24 


F15386+3807 


FIRST J154031. 0+375814 


0, 


,1828 


11, 


.62 


-1 


,00 


100 


.00 


ff 


1 


F15391+3214;SE 


2MASX J15410590+3204466 


0, 


,0530 


11, 


.39 


-24, 


,06 


100, 


.00 


wgs 


24 


F15394+3532 


F15394+3532 


0. 


.1235 


11, 


.81 


-1, 


.00 


100, 


.00 


wgs 


24 


F15519+3537 


2MASX J15534897+3528042 


0. 


.0842 


11, 


.66 


-25. 


.54 


100, 


.00 


wgs 


24 


F15543+3013 


FIRST J155623.3+300443 


0, 


.121 


11, 


.71 


-1, 


.00 


100, 


.00 


wgs 


2 


F15543+4158:NW 


2MASX J15560483+4149304 


0. 


.134 


11, 


.74 


-25. 


.55 


74, 


.96 


wgs 


2 


F15543+4158:SE 


FIRST J155605.3+414922 


0. 


.1342 


11, 


.26 


-24. 


.53 


25, 


.03 


wgs 


24 


F15549+4201 


2MASX J15563641+4152501 


0, 


.0348 


11, 


.16 


-24, 


.30 


100 


.00 


wgs 


24 


F16030+2040 


2MASX J16051287+2032326 


0, 


,015 


11, 


.02 


-23, 


,38 


100, 


.00 


rbgs 


5 


F16104+5235:NE 


2MASX J16114086+5227270 


0, 


,0295 


11, 


.51 


-24, 


,70 


100, 


.00 


rbgs 


54 


F16333+4630:E 


F16333+4630E 


0, 


,1908 


11, 


.36 


-23, 


,20 


7, 


.49 


Ijy 


1 


F16474+3430:N 


F16474+3430N 


0, 


,1126 


11 


.43 


-23, 


,82 


15, 


.28 


Ijy 


4 


F16504+0228 


2MASX J16525886+0224035 


0, 


,0243 


11, 


.86 


-26, 


,02 


100, 


.00 


rbgs 


5 


F16577+5900:SE 


2MASX J16583138+5856102 


0, 


,019 


11, 


.26 


-24, 


,82 


100, 


.00 


rbgs 


5 


F18131+6820 


2MASX J18125527+6821484 


0, 


,020 


11, 


.23 


-24, 


,83 


100, 


.00 


rbgs 


6 


F18425+6036 


2MASX J18431242+6039121 


0, 


,013 


11, 


.06 


-24, 


,45 


100, 


.00 


rbgs 


5 


F20550+1655:SE 


2MASXi J2057240+170735 


0, 


,035 


11, 


.88 


-23, 


,36 


100, 


.00 


rbgs 


3 


F22204-0214:NW 


2MASX J22225728-0159290 


0, 


,139 


11, 


.94 


-24, 


,30 


57, 


.27 


wgs 


2 


F22204-0214;SE 


... a- 


0, 


,140 


11, 


.82 


-23, 


,99 


42, 


.72 


wgs 


2 


F23488+2018 


2MASX J23511863+2034404 


0. 


,018 


11, 


.53 


-24. 


,58 


100, 


.00 


rbgs 


3 


ULIRGs 


F01298-0744 


F01298-0744 


0.1361 


12.41 


-24.38 


100.00 


Ijy 


1 


F03250+1606 


F03250+1606 


0.1290 


12.19 


-25.93 


100.00 


Ijy 


1 


F04313-1649 


F04313-1649 


0.2672 


12.72 


-25.14 


100.00 


Ijy 


1 
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TABLE 1 — Continued 



IRAS FSC 


Nuclear ID 


Rcdshift 


Lib. 


Ks 


% Lum. 


Sample 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 



F08136+3110 


FIRST J081645. 7+310119 


0, 


.4070 


12, 


.46 


-26, 


.30 


100.00 


ff 


1 


F08201+2801 


F08201+2801 


0, 


.1678 


12, 


.37 


-25. 


.28 


100.00 


ijy 


4 


F08208+3211 


FIRST J082354.6+320212 


0, 


.3955 


12, 


.49 


-25. 


.83 


100.00 


s 


1 


F08474+1813 


F08474+1813 


0, 


.1454 


12, 


.28 


-24, 


.18 


100.00 


Ijy 


1 


F08572+3915:NW 


IRAS 08572+3915NW 


0, 


.058 


12, 


.04 


-24, 


.09 


86.43 


ijy 


46 


F08591+5248 


F08591+5248 


0, 


.1573 


12, 


.30 


-26. 


.03 


100.00 


Ijy 


4 


F09039+0503 


F09039+0503 


0, 


.1252 


12, 


.16 


-25, 


.18 


100.00 


Ijy 


14 


F09116+0334:W 


F09116+0334W 


0, 


.1454 


12, 


.22 


-26. 


.35 


94.61 


Ijy 


14 


F09245+3300 


FIRST J092734.7+324713 


0, 


.222 


12, 


.30 


-1, 


.00 


100.00 


wgs 


2 


F09539+0857 


F09539+0857 


0, 


.1290 


12, 


.19 


-24, 


.37 


100.00 


Ijy 


1 


F10091+4704 


F10091+4704 


0, 


.2451 


12, 


.70 


-25. 


.96 


100.00 


Ijy 


1 


F10156+3705 


FIRST J101834.5+364951 


0, 


.4895 


12, 


.82 


-26. 


.08 


100.00 


s 


1 


F10494+4424 


F10494+4424 


0, 


.0919 


12, 


.23 


-25, 


.05 


100.00 


Ijy 


4 


F10565+2448 


2MASX J10591815+2432343 


0, 


.0430 


12, 


.04 


-25. 


.03 


100.00 


rbgs 


3 


F11028+3130 


F11028+3130 


0, 


.1986 


12, 


.49 


-24, 


.65 


100.00 


Ijy 


4 


F11387+4116 


F11387+4116 


0, 


.1489 


12, 


.26 


-25, 


.50 


100.00 


Ijy 


14 


F11506+1331 


F11506+1331 


0, 


.1274 


12, 


.39 


-25. 


.39 


100.00 


Ijy 


4 


F12112+0305:NE 


F12112+0305NE 


0, 


.073 


12, 


.00 


-24, 


.21 


45.64 


Ijy 


436 


F12112+0305:SW 


F12112+0305SW 


0, 


.0730 


12, 


.08 


-24, 


.40 


54.36 


Ijy 


43 


F13469+5833:W 


F13469+5833W 


0, 


.1575 


12, 


.04 


-24, 


.64 


59.11 


Ijy 


4 


F14348-1447:SW 


F14348-1447SW 


0, 


.083 


12, 


.15 


-25, 


.18 


62.18 


Ijy 


6 


F15206+3342 


F15206+3342 


0, 


.1254 


12, 


.26 


-25. 


.61 


100.00 


Ijy 


42 


F15250+3608 


2MASX J15265942+3558372 


0, 


.0552 


12, 


.02 


-24, 


.27 


100.00 


rbgs 


46 


F15414+3238 


FIRST J154326. 4+322859 


0, 


.2033 


12, 


.12 


-1, 


.00 


100.00 


wgs 


24 


F16007+3743 


2MASX J16023279+3734532 


0, 


.185 


12, 


.09 


-26. 


.18 


100.00 


wgs 


2 


F16333+4630:W 


F16333+4630W 


0, 


.1908 


12, 


.45 


-25, 


.85 


92.51 


Ijy 


1 


F17068+4027:E 


F17068+4027E 


0, 


.1794 


12, 


.34 


-25, 


.31 


78.04 


Ijy 


1 


F22381-1337 


2MASX J22404826-1321360 


0, 


.110 


12, 


.02 


-24, 


.64 


100.00 


wgs 


2 


F22491-1808 


F22491-1808 


0, 


.076 


12, 


.17 


-24, 


.92 


100.00 


Ijy 


26 


F23234+0946:W 


F23234+0946W 


0, 


.1279 


12, 


.07 


-25, 


.41 


78.72 


Ijy 


1 


F23365+3604 


2MASX J23390127+3621087 


0, 


.065 


12, 


.13 


-25. 


.22 


100.00 


rbgs 


2 



Ref eren ces. — 1 = IRupke et al.l ll200l. [2005bh : 2 = IKim et all l(T99l) : 3 = IWu et al.1 iflggl l: 4 = SDSS; 5 = IMoustakas k. KennicuttI 
1(2003] : 6 = ILiu fc Kennicutti 

Note. — Col.(l): Galaxy nucleus designation from the IRAS Fa int Source Catalog. Two sources are fflylS sources but not in the 
FSC, and a third is not an IRAS source (Fl_5; iClemen ts et al.H2001l). The nucleus is specified if the IRAS flux is unresolved or spread 
among two or more nuclei. Col. (2): Designation, in NED nomenclature, which precisely specifics in sky coordinates t he galaxy nucleu s 
from which the IR AS flux and spectrum originate. For 1 Jy ULIRGs, the nuclei are already known (Kim'et al. 2002; VeilleiD retaIll200g) . 
For non-1 Jy ULIRGs and all LIRGs, we specify a unique galaxy nucleus desig n ation. 2M ASS sources are used where possible. Col. (3): 
Redshift. Sources are, in order of preference, the SDSS, IRupke et aTl l l2005lj) . IKim et alj 819951) . iWu et al.l ai998i') . or NED. Col.(4): 
8 — 1000 fim infrared luminosity. The total galaxy lu mino sity has been multiplied by the factor in Column 5. Col. (5): A's-band absolute 
magnitude, compiled and computed as described in i|4.1l Col. (6); Fraction of (unresolved) total infrared luminosity originating in this 
nucleus, based on j-C-b and flux ratios for individual n uclei. Col. (7): Sample from which galaxy is taken: rb gs = Revised Bright Galaxy 
Sample (|Sanders et al.l'2003), liy = the 1 Jy sample IIKim fc Sanders! [T998D, wgs = Wa rm Galaxy Sample I IKim et al.nl995i'l . 2jy = the 
2 Jy sample (Strauss et al. 1993), ^^id ff = the FIRST /FSC sample llStanford et al.1120001 '). Col. (8): Spectrum reference; the first reference 
is the data we use in our analysis. 

^ A 2MASS designation exists for this galaxy, but it is attached to another nucleus. 
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TABLE 2 
New Line Fluxes 



IRAS FSC 


E(B-y) 


W|-(Ha) 




[Oil] 


H/3 


[OIII] 


[OI] 


Ha 


[N II] 


[SII] 


[SII] 






37264-3729 


4861 


5007 


6300 


6563 


6583 


6716 


6731 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 




LlKCis 


F01173-I-1405 


0.93 


-202.00 


3.29 


9.83e-15 


7.17e-15 


6.67e-15 


2.73e-15 


6.02e-14 


2.50e-14 


1.05e-14 


8.90e-15 


F01364-1042 


1.28 


-18.30 


9.84 


1.83e-15 


3.63e-16 


3.94e-16 


7.32e-16 


3.19e-15 


3.08e-15 


1.27e-15 


9.90e-16 


Fl_5 


0.91 


-119.48 


6.13 


1.49e-16 


1.43e-16 


5.56e-17 




1.14e-15 


5.13e-16 


1.86e-16 




F08572-|-3915:SE 


0.44 


-45.39 


4.02 


1.66e-15 


8.70e-16 


4.80e-16 


1.20e-16 


3.89e-15 


1.15e-15 


7.87e-16 


5.89e-16 


F09126+4432:NE 


1.45 


-14.40 


8.72 


1.83e-15 


6.29e-16 


3.62e-16 


3.80e-16 


6.53e-15 


6.14e-15 


1.47e-15 


1.07e-15 


F09209+3943:W 


1.07 


-11.03 


5.10 


8.19e-16 


2.00e-16 


2.09e-16 




1.45C-15 


9.87C-16 


4.85C-16 


3.52C-16 


F09320+6134 


1.63 


-35.72 


9.80 


1.95e-15 


1.02e-15 


1.67C-15 


1.07C-15 


1.60C-14 


1.92C-14 


3.34C-15 


2.53C-15 


F10190+1322:W 


0.89 


-36.10 


6.14 


9.21e-16 


1.31e-15 


2.16C-16 


2.35C-16 


9.31C-15 


5.47C-15 


1.08C-15 


9.88C-16 


F10190+1322:E 


1.88 


-18.90 


10.17 


5.05e-16 


1.57e-16 


1.86C-16 


3.43C-16 


2.80e-15 


2.41C-15 


6.55C-16 


4.80e-16 


F10594+3818;SW 


0.51 


-98.90 


9.76 


4.480-15 


3.10C-15 


1.60C-15 


6.97C-16 


1.49e-14 


7.74G-15 


2.62e-15 


2.05e-15 


F12592+0436 


1.34 


-17.02 


9.03 


1.25e-15 


4.46c- 16 


3.81C-16 


5.46e-16 


4.24e-15 


3.63C-15 






F13136+6223 


0.63 


-244.90 


4.24 


5.57e-14 


3.00e-14 


3.51e-14 


5.64e-15 


1.78C-13 


6.93C-14 


2.27C-14 


2.10C-14 


F15043+57ryl:S 


0.52 


-80.16 


5.82 


2.36C-15 


1.27C-15 


7.97c- 16 


2.58C-16 


6. 28c- 15 


8.18C-16 


1.18C-15 


9.12C-16 


F15386+3807 


1.21 


-75.33 


4.84 


1.56C-16 


9.91C-17 


8.59C-17 


4. 70c- 17 


l.llc-15 


7. 54c- 16 


2.00C-16 


1.56C-16 


F16333+4630;E 


0.49 


-54.53 


5.42 


2.10e-16 


9.50e-17 


8.71e-17 


2.50e-17 


4.47e-16 


1.73e-16 


8.70e-17 


6.31e-17 


F16474+3430:N 


0.57 


-66.99 


7.80 


2.12e-15 


1.02e-15 


6.19e-16 


1.70e-16 


5.20e-15 


7.02e-16 


9.81e-16 


7.10e-16 



ULIRGs 



F01298-0744 


0.92 


-39.24 


3.60 


4.14e-16 


2.26e-16 


2.33e-16 


1.80c- 


16 


1.76e-15 


9.46C-16 






F03250+1606 


0.90 


-23.10 


6.95 


7.92e-16 


5.54e-16 


5.65C-16 


5.69c- 


16 


3.61e-15 


3. 64c- 15 






F04313-1649 


1.66 


-38.68 


8.67 


3.84C-16 


2.23C-16 


1.85c- 16 


1.22c- 


16 


3.75C-15 


9.30C-16 


3.54C-16 


2.71C-16 


F08136-I-3110 


0.75 


-118.41 


6.72 


5.90C-16 


1.88C-16 


2.10e-16 


1.35e- 


16 


1.23c- 15 


5.92C-16 


3.40C-16 


2.00C-16 


F08201+2801 


0.81 


-80.65 


7.71 


1.88C-15 


8.30C-16 


6.51C-16 


3.07c- 


16 


5.71C-15 


2. 49c- 15 


1.06C-15 


8.60C-16 


F08208+3211 


0.79 


-95.11 


2.69 


2.67C-16 


1.25C-16 


1.41C-16 


4.80c- 


17 


8.83C-16 


4. 74c- 16 




1.42C-16 


F08474-I-1813 


1.56 


-11.77 


6.75 


1.57e-16 


3.30e-17 


6.39e-17 


5.15e- 


17 


3.94e-16 


5.44e-16 


1.30e-16 


1.09e-16 


F08572-|-3915:NW 


0.69 


-59.57 


7.00 


2.09e-15 


1.25e-15 


1.12e-15 


4.60e- 


16 


7.29e-15 


2.89e-15 


1.85e-15 


1.35e-15 


F0859H-5248 


0.82 


-42.95 


9.49 


1.96e-15 


9.42e-16 


4.55e-16 


4.32e- 


16 


5.91e-15 


4.40e-15 


1.22e-15 


9.98e-16 


F09039-I-0503 


0.54 


-85.37 


8.88 


2.75e-15 


2.07e-15 


1.15e-15 


1.03e- 


15 


1.03e-14 


7.52e-15 


2.68e-15 


2.56e-15 


F09116-|-0334:W 


1.06 


-25.11 


11.69 


2.30e-15 


l.OOe-15 


5.46e-16 


6.82e- 


16 


7.15e-15 


7.50e-15 




1.35e-15 


F09539-I-0857 


0.97 


-33.12 


9.30 


1.82e-15 


3.19e-16 


6.69e-16 


4.88e- 


16 


2.30e-15 


3.65e-15 


9.83e-16 




F10091-I-4704 


1.09 


-25.28 


5.86 


2.87e-16 


6.95e-17 


8.72e-17 


1.06e- 


16 


6.00e-16 


6.24e-16 


2.68e-16 


1.94e-16 


F10156-I-3705 


0.21 


-229.17 


8.46 


7.44C-16 


2.56e-16 


3.75C-16 






9.08e-16 


3.50e-16 




1.30e-16 


F10494+4424 


1.43 


-50.70 


7.00 


8.37C-16 


3.06c- 16 


4.72C-16 


4.62e- 


16 


4.21C-15 


1.07C-15 


1.07C-15 


8.92C-16 


F11028+3130 


-0.01 


-8.44 


9.09 


4.01C-16 


2.82C-16 


8.91C-17 






4.17e-16 


3.77e-16 


1.90e-16 


1.05C-16 


F11387+4116 


1.02 


-27.00 


4.34 


1.13e-15 


4.95C-16 


3. 10c- 16 


4.00c- 


16 


4.13e-15 


4.10e-15 


1.16e-15 


l.OOc-15 


F11506+1331 


1.13 


-78.60 


2.99 


1.77C-15 


6. 16c- 16 


8.80c- 16 


5.76c- 


16 


6.41e-15 


2.74e-15 


1.43e-15 


1.13C-15 


F12112+0305:NE 


0.89 


-76.70 


8.86 


4.63C-15 


1.48C-15 


2.61C-15 


1.07c- 


15 


1.09C-14 


4.62e-15 


2.74e-15 


1.99e-15 


F12112+0305:SW 


1.33 


-38.30 


9.17 


1.15C-15 


3.55C-16 


3.82C-16 


4.93c- 


16 


4.00C-15 


2.60C-15 


1.09e-15 


8.07e-16 


F13469+5833:W 


0.98 


-26.02 


6.44 


6.73C-16 


3.01C-16 


2.17C-16 


1.25c- 


16 


2. 25c- 15 


1.86C-15 


4.80e-16 


3.44e-16 


F15206+3342 


0.58 


-401.00 


0.00 


3.01C-14 


1.93C-14 


5.41C-14 


3.38c- 


15 


1.09C-13 


2.04C-14 


9.54e-15 


7.87C-15 


F15250+3608 


0.75 


-47.65 


5.44 


7.03C-15 


3.11C-15 


3.87C-15 


1.49c- 


15 


1.95G-14 


9.34e-15 


4.13C-15 


3.32C-15 


F16333-|-4630:W 


0.82 


-113.20 


6.31 


1.13e-15 


5.52e-16 


8.07e-16 


3.58e- 


16 


3.92e-15 


2.46e-15 


8.72e-16 


7.13e-16 


F17068-|-4027;E 


0.35 


-187.26 


11.34 


1.25e-15 


8.40e-16 


1.22e-15 


2.21e- 


16 


3.44e-15 


1.07e-15 


6.70e-16 


5.40e-16 


F23234-|-0946:W 


0.65 


-167.52 


0.00 


4.00e-15 


1.45e-15 


1.68e-15 


1.02e- 


15 


8.90e-15 


9.90e-15 


2.66e-15 





Note. — - Col.(l): Galaxy designation from the IRAS Faint Source Catalog. Col. (2): Optical extinction, measured from the Balmer decrement. 
Col. (3): Ha emission line equivalent width, in A, uncorrected for underlying absorption. Col. (4): H/3 stellar absorption line equivalent width. 
Col. (5-12): Observed line fluxes, in erg cm~^. Fluxes are uncorrected for extinction or underlying stellar absorption, except for H/3 (which 
has had stellar absorption contamination removed). Line flux measurement errors range are 5% for strong lines, and up to 20 — 30% in a few 
individual cases for weak lines (e.g., [O I], [O III], and/or [S II]). 
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TABLE 3 

Abundances 



IRAS FSC 


log(R23) 


log(032) 






12-log(0/H) 












1 


2 


3 


4 


5 


6 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




LIRGs 


F00189+3748 


0.49 


-1.72 


8.09 


8.82 


8.93 


9.09 


8.95 


9.14 


F00267+3016:NW 


1.36 


-1.60 














F00267+3016:SE 


0.91 


-1.01 


7.70 


8.33 


8.44 


8.68 


8.84 


9.10 


F01173+1405 


0.71 


-0.92 


8.07 


8.63 


8.71 


8.75 


8.88 


8.92 


F01364— 1042 


1.37 


-1.84 














F01484+2220 


0.33 


-1.10 


8.38 


8.94 


9.04 


8.98 


9.04 




F02248+2621 


0.34 


-0.86 


8.45 


8.94 


9.03 


8.91 


9.13 


9.40 


F02512+1446 


0.65 


-1.14 


8.07 


8.69 


8.79 


8.85 


8.92 


9.07 


Fl_5 


0.54 


-1.25 


8.16 


8.79 


8.89 


8.92 


8.95 


9.06 


F04315-0840 


0.32 


-0.15 


8.69 


8.96 


9.05 


8.78 


9.17 


9.31 


F06538+4628:SW 


0.71 


-1.24 


7.95 


8.61 


8.72 


8.86 


8.88 




IRAS07062+2041 


0.33 


-0.78 


8.49 


8.95 


9.04 


8.90 


9.06 


8.97 


IRAS07063+2043 


0.72 


-1.41 


7.88 


8.58 


8.70 


8.89 


8.94 




F07256+3355 


0.94 


-2.06 


7.32 


8.13 


8.39 


8.99 


8.75 


9.27 


F08572+3915:SE 


0.59 


-1.05 


8.17 


8.75 


8.84 


8.85 


8.84 


8.76 


F09046+1838 


0.76 


-1.43 


7.81 


8.52 


8.65 


8.87 


8.84 


9.12 


F09126+4432:NE 


1.21 


-2.02 














F09209+3943:E 


0.48 


-1.27 


8.21 


8.84 


8.94 


8.96 


8.99 


9.09 


F09209+3943:W 


1.19 


-1.59 














F09218+3428 


0.58 


-0.90 


8.24 


8.77 


8.85 


8.81 


8.95 


9.01 


F09320+6134 


1.16 


-1.29 














F09333+4841 


0.90 


-1.24 


7.62 


8.32 


8.45 


8.76 


8.84 


9.13 


F09338+3133 


0.43 


-1.59 


8.18 


8.87 


8.98 


9.09 


8.99 


9.19 


F09339+2835 


0.56 


-1.09 


8.19 


8.78 


8.87 


8.87 


8.96 


9.05 


F09399+2830 


0.97 


-0.94 


7.60 


8.23 


8.34 


8.64 


8.77 


8.97 


F10190+1322:W 


0.34 


-1.51 


8.27 


8.94 


9.04 


9.10 


9.06 


9.36 


F10190+1322:E 


1.47 


-1.97 














F10203+5235 


0.71 


-0.50 


8.26 


8.66 


8.72 


8.66 


9.01 


8.99 


F10594+3818:SW 


0.51 


-0.98 


8.28 


8.83 


8.91 


8.87 


9.00 


9.09 


F11231+1456 


0.52 


-1.77 


8.05 


8.79 


8.90 


9.09 


8.99 


9.32 


F12116+5448 


0.55 


-0.36 


8.48 


8.82 


8.88 


8.71 


9.04 


9.15 


F12251+4026 


0.59 


-1.29 


8.09 


8.74 


8.84 


8.92 


8.97 




F12592+0436 


1.15 


-1.69 














F13136+6223 


0.73 


-0.74 


8.13 


8.62 


8.69 


8.70 


8.87 


9.03 


F13182+3424 


1.16 


-1.52 














F13373+0105:NW 


0.02 


0.17 


8.84 


9.06 


9.18 


8.79 


9.23 


9.45 


F13373+0105:SE 


0.58 


-0.85 


8.26 


8.77 


8.85 


8.81 


8.67 




F13458+1540:NE 


0.68 


-1.23 


8.00 


8.65 


8.75 


8.86 


8.85 


8.96 


F13458+1540:SW 


0.44 


-0.68 


8.45 


8.89 


8.97 


8.83 


9.04 


9.04 


F14179+4927 


0.31 


-1.27 


8.35 


8.96 


9.05 


9.03 


9.13 




F14280+3126 


0.37 


-1.49 


8.24 


8.92 


9.02 


9.08 


8.99 


9.29 


F14416+6618 


0.77 


-0.30 


8.27 


8.61 


8.65 


8.59 


8.87 


8.85 


F15043+5754:S 


0.62 


-1.02 


8.15 


8.72 


8.81 


8.83 


8.80 


8.36 


F15163+4255:NW 


0.69 


-0.87 


8.12 


8.66 


8.74 


8.75 


8.93 




F15359+3139 


0.90 


-1.60 


7.51 


8.28 


8.46 


8.86 


8.83 


9.28 


F15384+3841 


0.51 


-1.09 


8.24 


8.83 


8.92 


8.90 


8.94 


9.37 


F15386+3807 


0.87 


-1.24 


7.68 


8.37 


8.50 


8.76 


8.88 


9.18 


F15391+3214:SE 


0.71 


-1.78 


7.81 


8.57 


8.72 


9.01 


8.88 


9.28 


F15394+3532 


0.49 


-0.84 


8.35 


8.85 


8.93 


8.84 


8.94 


8.95 


F15519+3537 


0.85 


-1.30 


7.70 


8.41 


8.54 


8.79 


8.84 


9.27 


F15543+3013 


0.59 


-0.89 


8.23 


8.76 


8.84 


8.80 


9.00 


9.17 


F15543+4158:NW 


0.57 


-0.69 


8.34 


8.79 


8.87 


8.77 


8.95 


8.72 


F15543+4158:SE 


0.68 


-0.04 


8.45 


8.72 


8.74 


8.59 


8.87 


8.84 


F15549+4201 


0.76 


-0.69 


8.10 


8.58 


8.65 


8.68 


8.91 


9.04 


F16030+2040 


0.64 


-0.78 


8.22 


8.72 


8.79 


8.76 


8.80 


8.68 


F16104+5235:NE 


0.38 


-0.82 


8.44 


8.92 


9.01 


8.89 


9.06 




F16333+4630:E 


0.71 


-0.88 


8.10 


8.64 


8.72 


8.75 


8.86 


8.90 


F16474+3430:N 


0.68 


-1.15 














F16504+0228 


1.05 


-1.28 














F16577+5900:SE 


0.61 


-1.29 


8.07 


8.73 


8.83 


8.91 


8.95 


9.07 


F18131+6820 


0.23 


-0.29 


8.69 


9.00 


9.09 


8.85 


9.15 


9.38 


F18425+6036 


0.45 


-0.97 


8.33 


8.87 


8.96 


8.89 


9.10 


9.23 


F20550+1655:SE 


(1.8.1 


-0.31 


8.20 


8.55 


8.58 


8.58 


8.74 


8.90 


F22204-0214:NW 


0.98 


-1.42 


7.36 


8.13 


8.31 


8.77 


8.76 


9.12 


F22204-0214:SE 


0.42 


-1.31 


8.25 


8.89 


8.98 


9.00 


9.01 


9.40 


F23488+2018 


0.60 


-1.26 


8.09 


8.74 


8.84 


8.90 


8.97 


9.23 


ULIRGs 


F01298-0744 


0.82 


-1.02 


7.87 


8.48 


8.58 


8.73 


8.86 




F03250+1606 


0.72 


-0.86 


8.08 


8.62 


8.70 


8.73 


9.02 
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1 


2 


3 


4 


5 


6 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


r U4oiv5- 


1 «A O 

— iD4y 


111 
1.11 


1 «^ 
-i.DO 














r UoioD- 


i~o±iu 


u.yo 


1 1 R 
- 1 . ID 


7 

1 . 00 


8 9A 
o. Z^ 


o.Oo 


8 71 
o. / 1 


8 7« 
o. / D 


8 88 

O.OO 


X" UoZUi- 


1 osni 

i~Zf5U± 


U.oo 


1 99 
- 1 . ZZ 


7 7fi 
1 . / D 


8 AA 


O.OD 


8 78 
o. / o 


8 70 

o. / y 


8 Q7 

o.y / 


r UoZUo- 


-j-ozi ± 


u.oo 


-u.yD 


7 87 

1 .O 1 


8 A7 
0.4 / 


O.OD 


8 71 
o. / 1 


8 87 

O.O 1 


8 QQ 

o.yy 


r Uo't ( '4:- 


"h iolo 


1 AQ 


- 1 .Do 














r UoO / Z- 


"hoy 10. IN vv 


n RS 

U.Do 


n 87 

-U.O 1 


o. 10 


O.DD 


8 7A 
O. 1 4 


O. / D 


8 87 

O.O 1 


8 7^ 
o. 1 


r UoOyi' 


1 croAS 

i-oz^o 


u. / O 


1 Afi 
- 1 .^D 


7 78 
(.to 


8 

O. OU 


O.DO 


8 88 

O.OO 


8 QI 

o.y 1 


Q 99 

y .ZZ 


FTlQfl'iQ 

r uyuoy- 


i~UOUo 


U.Ol 


n Qi 
-u.yi 


8 "^l 
o.Ol 


8 8*^ 
o.oO 


8 Q9 

o.yz 


8 9.^ 
O.oO 


Q 07 

y.u 1 


Q 01 

y .ui 


r uy i ID" 


TKJOO^. vv 


u.yo 


- 1 .DO 


7 A*^ 
1 .40 


8 91 
o. Z 1 


8 AO 
0.4U 


8 8fl 
O.oD 


8 8fl 
O.oD 


Q "^1 

y .Ol 


r uyz^O" 


i-oouu 


U.DZ 


-U.oD 


8 91 
o. Zl 


8 7*^ 
O. /O 


8 81 
O.Ol 


8 77 

O. / 1 


8 QQ 

o.yy 


Q 1 

y . lo 


-T uyooy- 


i-UoO / 


l.OU 


-1 .ou 














r luuyi- 


"1-4 i U4 


1 91 
l.Zl 


1 ^^ 

- 1 .Ol 














■pi ni pifi 




n 7"^ 
u. / o 


-U. OD 


8 18 
o. lo 


O.Dl 


O.DD 


O.DD 


8 88 

O.OO 


8 88 

O.OO 


■pi OAQA 

X" lu^y^- 


1/1/10/1 


1 91 
l.Zl 


1 AO 














■pi nf^fi'^ 


i_0/1 AS 


n Qi 
u.y 1 


1 9A 

- 1 . Z'y: 


7 fin 

1 .DU 


8 "^l 

O.Ol 


8 AA 
0.4^ 


8 7^ 
o. / 


8 8(^ 
O.oD 


Q 99 

y .ZZ 


■pi 1 098 


"hOlOU 


n 9fi 

U.ZD 


n 87 

-U.O 1 


8 f^n 
o. ou 


8 Q8 

o.yo 


Q 08 

y.uo 


8 07 

o.y 1 


Oil 

y. 1 1 


8 QQ 

o.yy 


■pi 1 QQ7 

r i loo 1 - 


J_A1 1 R 
■f-^l ID 


n Q9 
u.yz 


1 ^^ 

- 1 .Ol 


7 f^n 
/ . ou 


8 9fi 
O. ZD 


8 A'X 
0.40 


8 8*^ 
O.oO 


8 QI 

o.y 1 


Q 1 8 

y . lo 


pi 1 c;fiR 


lool 


1 no 
i.uy 


1 9A 














pi Olio 




1 OA 


n QQ 
-u.yy 














■pi Olio 


i-uouo-o vv 


1 91 
l.Zl 


- 1 .DO 














■pi Q/lflQ 


i-oooo- VV 


n Qn 
u.yu 


1 'iQ 

- 1 .oy 


7 ^^8 
/ . Oo 


8 "^1 
O.Ol 


8 A(\ 
0.4D 


8 80 
o.oU 


8 87 

O.O 1 


y .ZD 


■pi /I Q/IQ 


— i.^^ i .O VV 


n f^7 

U.O ( 


U.Ul 


8 ^7 
o.Ol 


8 89 
o.oZ 


8 8(^ 
O.oD 


8 (\A 
O.D4 


Q 1 9 

y. iz 


8 Qf^ 

o.yo 


F15206- 


f3342 


0.82 


-0.09 


8.28 


8.56 


8.57 


8.44 


8.76 


8.80 


F15250- 


f3608 


0.84 


-0.90 


7.88 


8.45 


8.54 


8.69 


8.83 


8.92 


F15414- 


f3238 


0.73 


-1.21 


7.93 


8.58 


8.69 


8.83 


8.92 




F16007- 


f3743 


1.06 


-1.46 














F16333- 


f4630:W 


0.85 


-0.80 


7.91 


8.45 


8.53 


8.66 


8.90 


9.01 


F17068- 


f4027:E 


0.61 


-0.27 


8.45 


8.77 


8.82 


8.66 


8.94 


8.70 


F22381- 


-1337 


0.61 


0.01 


8.53 


8.78 


8.82 


8.52 


8.95 


8.82 


F22491- 


-1808 


0.57 


-0.87 


8.26 


8.78 


8.86 


8.81 


8.96 


9.08 


F23234- 


f0946:W 


0.87 


-0.99 


7.79 


8.41 


8.51 


8.71 


9.01 


9.21 


F23365- 


f3604 


0.64 


-1.14 


8.08 


8.70 


8.79 


8.85 


9.00 


9.03 



Note. — Col.(l): Galaxy designation from the IRAS Faint Source Catalog. Col.(2): Logarithm of R23 = {/([O II]AA3726. 3729) + 
/([O III]AA4959, 5007)}//(///3). Col.(2): Logarithm of O32 = f(\0 III]AA4959, 5007)/ /([O II1AA3726, 3729). Col.(4-9): Oxygen 
abun danc e from six different d iagnostic/ca libration systems. 1 = IPilvugin Thuanl 120051 : 2 =jMcGaugh 1991; Kuzio dc Narav ct aLl 
[200I : 3 = ITremonti et al.ll200l : 4 and 6 = ICharlot fc LonghettH 120011 . Cases F and A: and 5 = lKewlev fc Dopitali2002i . See ^3l1 for a 
description of each system. 



